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HICHLICHT  SUt-II'lARY 

A  major  detrimental  aspect  of  the  traditional  urban  transportation 
planning  modelling  process,  especially  for  the  smaller  urban  areas,  has 
been  its  voracity  for  large  amounts  of  data.   The  collection  and  pro- 
cessing of  data  for  urban  travel  demand  modelling  is  both  costly  and 
time  consuming,  resulting  in  the  final  transportation  plan  being  com- 
pleted years  after  the  initial  data  collection  and  often  without  adequate 
opportunity  for  alternative  testing.   In  order  to  possibly  remedy  this 
situation  of  high  costs  and  slow  turnaround  time,  various  synthetic 
modelling  techniques  have  been  suggested  and  tested  in  recent  years 
with  mixed  results.   As  a  part  of  a  detailed  study  of  synthetic  travel 
demand  modelling  techniques  that  can  be  applied  to  small  urban  areas 
(population  50,000  to  250,000)  in  Indiana,  the  present  study  concentrated 
solely  on  the  trip  distribution  phase  of  the  modelling  process. 

The  most  commonly  used  trip  distribution  model  is  the  gravity  model 
in  which  calibration  of  the  model  primarily  involves  the  establishment 
of  the  appropriate  travel  time  impedance  function.   The  present  study 
evaluated  the  performance  of  some  of  the  various  suggested  synthetic 
gravity  model  calibration  techniques  and  developed  an  improved  procedure 
for  synthetic  self-calibration  of  the  gravity  model.   The  primary  purpose 
of  this  study  was  to  eliminate  the  need  for  the  home-interview  survey 
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from  which  calibration  data  is  normally  obtained  and  to  do  so  without 
sacrificing  accuracy  of  the  results. 

Data  from  a  wide  range  of  small  urban  areas  were  used  in  deter- 
mining variability  of  different  characteristics  among  the  small  urban 
areas  in  order  to  evaluate  the  possible  application  of  borrowed  travel 
time  impedance  factors,  or  their  associated  mathematical  functions, 
in  urban  areas  with  similar  characteristics.   Additionally,  the  com- 
plete travel  demand  modelling  data  from  three  small  urban  areas  of 
Indiana  were  used  as  test  cases  to  which  various  synthetic  trip  distri- 
bution modelling  techniques  were  applied  and  evaluated.   The  data  from 
these  three  urban  areas  also  provided  a  data  source  for  the  develop- 
ment of  the  recommended  synthetic  self-calibrating  trip  distribution 
model. 

The  results  of  this  research  effort  revealed  the  great  variabil- 
ity among  small  urban  areas  as  well  as  their  individual  uniqueness. 
Due  to  this  variability  and  uniqueness,  it  was  concluded  that  the 
most  commonly  suggested  and  used  synthetic  trip  distribution  modelling 
technique  of  borrowing  travel  time  impedance  values  from  similar  urban 
areas  cannot  provide  consistent  and  reliable  results.   The  same  con- 
clusion was  also  reached  for  the  technique  of  using  standardized 
travel  time  impedance  values.   In  addition,  the  present  study  also 
indicated  that  small  urban  areas  may  differ  greatly  in  terms  of  their 
sensitivity  to  travel  time  impedance  value  errors  when  measured  at  the 
traffic  assignment  level.   That  is,  one  small  urban  area's  models  may 
produce  acceptable  results  within  a  wide  travel  time  impedance  value 
range  while  another  small  urban  area's  models  may  provide  acceptable 
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results  only  within  a  narrow  band  of  travel  time  impedance  values. 
These  findings  indicated  the  need  for  an  improved  synthetic  trip  dis- 
tribution procedure  that  would  provide  consistent  and  reliable  results. 

The  development  of  the  recommended  synthetic  self-calibrating 
gravity  model  is  based  primarily  on  the  origin  zone  opportunity  mean 
travel  time,  that  is,  the  origin  zone  mean  travel  time  that  would 
occur  if  the  travel  time  impedance  values  were  constant.   Origin  zone 
opportunity  mean  travel  times  are  themselves  closely  related  to  the 
spatial  location  of  the  zones  and  the  land  use  activities  within  the 
area.   These  origin  zone  opportunity  mean  travel  times  are  then  used 
to  estimate  the  actual  origin  zone  mean  travel  times  and  the  origin 
zone  travel  time  impedance  function  parameter.   It  is  felt  that  such 
a  model  is  more  responsive  to  the  location  of  activities  than  the 
traditionally  calibrated  gravity  model. 

Application  of  the  recommended  model  to  three  different  urban 
areas  proved  successful  in  generating  acceptable  traffic  assignments 
when  compared  to  the  Origin-Destination  survey  assignment.   For  one 
test  case,  total  volume  group  percent  root-mean-square-error  was 
actually  reduced  from  67.8  percent  to  47.7  percent  when  compared  to 
the  traditionally  calibrated  gravity  model  traffic  assignment.   The 
recommended  model  thus  appears  to  be  capable  of  significantly 
decreasing  the  time  and  effort  presently  expended  in  trip  distribution 
modelling. 


CHAPTER  I 
INTRODUCTION 

Urban  areas  approaching  the  50,000  population  level  face  the 
formidable  task  of  initiating  a  comprehensive  transportation  planning 
process.   Other  urban  areas,  having  already  completed  the  initial  major 
transportation  planning  process,  need  to  continuously  update  plans  as 
the  urban  environment  and  the  needs  of  the  community  change.   Both  the 
initiating  and  the  updating  of  a  comprehensive  plan  require  considerable 
resources  of  time  and  money.   Furthermore,  the  emphasis  in  urban  trans- 
portation planning  is  currently  changing  from  traditional  long-range 
system  planning  to  short-range,  quick  response  improvement  programs 
with  the  aim  of  optimizing  the  utility  of  the  existing  facilities  (7). 
The  traditional  travel  forecasting  process  must  be  modified  if  it  is  to 
remain  a  relevant  and  coordinated  phase  of  both  the  short-range  and 
long-range  planning  process. 

Consequently,  it  is  necessary  to  make  efforts  toward  simplifying 
planning  techniques  such  that  savings  can  be  achieved  in  both  the  costs 
and  time  involved  in  planning.   Additionally,  simplified  planning 
techniques  would  allow  greater  attention  to  be  given  to  the  more 
immediate  needs  of  transportation  system  management  planning,  where  it 
is  imperative  that  quick  response  time  be  inherent  to  the  techniques  (16) 

In  this  regard,  a  recent  National  Cooperative  Highway  Research 
Project  (NCHRP)  study  (11)  investigated  the  ability  of  various  existing 


travel  estimation  procedures  to  provide  information  within  a  quick. 

response  framework  for  analyzing  urban  transportation  policy  issues. 

The  following  statement  in  regards  to  the  results  of  this  investigation 

was  made  in  the  summary  of  the  report: 

No  existing  travel  estimation  technique  was  found  adequate 
to  respond  to  all  the  many  questions  raised.   Because  of  the 
widespread  use  and  understanding  of  the  conventional  four-step 
(trip  generation,  distribution,  mode  split,  and  traffic 
assignment)  travel  simulation  procedure,  it  was  recommended 
that  this  procedure  continue  to  form  at  least  a  partial  basis 
on  which  the  transportation  planning  process  responds  to  new 
needs  in  the  future. 

The  author  of  the  report  further  recommended  that  a  set  of  capabilities 

be  developed  for  simplified  computerized  methods  for  the  conventional 

four-step  procedure  in  order  to  provide  efficient  and  quick  response 

to  transportation  planning  issues  at  the  regional  and  subarea  level. 

Purdue  University  and  the  Planning  Division  of  the  Indiana  State 
Highway  Commission  (ISHC)  recognized  the  need  for  simplified  procedures, 
and  therefore,  in  early  1977,  organized  a  research  project  to  evaluate 
the  use  of  synthetic  demand  modelling  techniques  for  transportation 
planning  in  the  small  urban  areas  of  Indiana.   The  Planning  Division  of 
ISHC  was  very  aware  of  the  past  resources  expended  in  the  transporta- 
tion planning  process  for  many  of  the  urbanized  areas  in  Indiana  and 
therefore  realized  that  the  results  derived  from  such  a  project  could 
be  very  beneficial. 

The  scope  of  the  overall  research  project  included  the  evaluation 
of  the  three  transportation  modelling  steps  of  trip  generation,  trip 
distribution,  and  traffic  assignment.   Because  the  results  of  this 
research  were  to  be  applicable  to  small  urban  areas  of  population  size 


50,000  to  around  250,000  persons,  it  was  felt  that  transit  planning 
could  be  modelled  separately,  thereby  eliminating  mode  split  from  the 
study.   This  same  conclusion  was  expressed  by  Grecco,  et .  al.  (6)  in  a 
study  conducted  for  the  NCHRP. 

In  the  first  phase  of  the  research  project,  conducted  by  Mahmassani 
(10),  synthetic  trip  generation  techniques  as  applicable  to  the  small 
urban  areas  of  Indiana  were  evaluated.   Mahmassani' s  research  indicated 
that  vehicle  trip  rates  per  household,  as  used  in  the  cross-classifica- 
tion technique,  were  sufficiently  stable  to  be  applicable  to  any  of  the 
urbanized  areas  within  Indiana.   The  cross-classification  technique 
is  well  documented  in  the  Federal  Highway  Administration's  (FHWA) 
publication  Trip  Generation  Analysis  (17). 

Origin-Destination  (0-D)  surveys  of  the  past  had  been  the  corner- 
stones from  which  trip  generation  models  were  originally  developed. 
The  transferrable  trip  rates  per  household  have  eliminated  the  need 
for  the  0-D  survey  for  trip  generation.   However,  0-D  surveys  had 
provided  other  information  critical  to  the  transportation  modelling 
process . 

Travel  time  frequency  distributions  developed  from  0-D  survey  trip 
tables  and  highway  network  skim  trees  (travel  time  matrix)  were  used  in 
the  calibration  of  the  trip  distribution  models  of  the  past.   Elimina- 
tion of  the  0-D  survey  would  remove  the  information  used  in  trip  distri- 
bution calibration.   Other  means  of  calibrating  the  model  would  have  to 
be  used  in  the  procedure. 


Trip  distribution,  using  synthetic  methods,  is  still  in  its  early 
development  and  testing  stage.   The  second  phase  of  the  synthetic 
modelling  project  therefore  focusses  on  trip  distribution,  specifically 
for  small  urban  areas  (population  size  range  from  50,000  to  about 
250,000  persons).   The  overall  objective  of  the  present  study  was  to 
reduce  the  quantity  of  time,  money,  and  manpower  consumed  in  trip 
distribution  modelling  without  adversely  affecting  the  accuracy  of 
the  results. 

This  research  project  examined  those  small  urban  areas  of  Indiana 
where  full  scale  transportation  studies  have  already  been  conducted,  in 
an  effort  to  make  more  efficient  use  of  the  information  already  made 
available  by  these  studies.  Data  from  outside  of  Indiana  was  also  used 
so  as  to  increase  the  sample  size.  In  addition,  the  complete  transpor- 
tation modelling  data  sets  from  three  urban  areas,  Lafayette,  Anderson, 
and  Muncie,  were  used  for  detailed  analysis  and  testing. 

The  work  plan  for  this  study  consisted  of  five  major  tasks: 

1.  Summarize  the  travel  time  statistics  of  the  Indiana  transpor- 
tation study  areas,  and  if  possible,  relate  these  statistics 
to  other  urban  area  characteristics.   Determine  if  these 
statistics  will  be  of  any  use  in  other  urban  areas  in  which 
the  transportation  planning  process  is  to  be  implemented  or 
updated. 

2,  Investigate  the  sensitivity  of  the  gravity  model  with  respect 
to  friction  factor  errors  introduced  by  the  use  of  borrowed 
or  standardized  friction  factors. 


3.  Investigate  the  possibility  of  using  areawide  travel  time 
statistics  and  other  variables  in  a  self-calibrating  gravity 
model. 

4.  Investigate  the  possibility  of  using  origin  zone  specific 
travel  time  statistics  in  a  synthetically  self-calibrating 
gravity  model. 

5.  Evaluate  and  make  recommendations  as  to  what  method  of 
synthetic  trip  distribution  modelling  should  be  used  in  the 
transportation  planning  process  for  small  urban  areas. 


CHAPTER  II 
TRAVEL  TIME  CHARACTERISTICS 

An  urban  area's  travel  time  characteristics  reflect  the  nature  of 
the  activity  system  of  the  urban  area  as  represented  by  the  land  use 
spatial  distribution  and  the  socio-economic  background  of  the  area. 
On  the  basis  of  a  given  activity  system  and  an  associated  transportation 
system,  decisions  are  made  by  individual  trip  makers  as  to  how  often, 
at  what  time,  to  what  place,  by  what  mode,  and  by  which  route  a  trip  is 
to  be  made.   The  resulting  cumulative  travel  pattern  is  generally 
quantified  by  trip  purpose  in  terms  of  the  mean  travel  time,  travel 
time  variance,  and  by  the  corresponding  travel  time  frequency  distribu- 
tion. 

Literature  Review 

Relationships  between  travel  time  characteristics  and  the  size  and 
structure  of  urban  areas  have  been  the  subject  of  several  studies  with 
the  general  conclusion  that  mean  travel  times  are  positively  related  to 
population  and  the  developed  land  area.   The  more  prominent  research 
efforts  in  this  subject  were  performed  by  Chan  (2)  and  Voorhees  et.al.(19) 

Voorhees,  et.  al.  extensively  studied  travel  time  (trip  length) 
characteristics  for  work  trips,  shopping  trips,  socio-recreational 
trips,  non-home  based  trips,  and  truck  trips  from  origin-destination 
data  collected  from  34  transportation  study  areas.   Study  area 


population  ranged  from  a  minimum  of  33,000  to  a  maximum  of  6,489,000 
persons. 

R.egression  analysis  was  used  by  Voorhees  et.  al.  to  test  the 

hypothesis  that  trip  length  (trip  duration)  increases  with  increasing 

population  size.   Equations  1  and  2  show  the  results  of  this  analysis 

for  the  home  based  work  (HBW)  and  non-home  based  (NHB)  trips, 
respectively. 


0.19 
Home  Based  Work  t  =  0.98  POP 

R^  =  0.71 


Non-Home  Based  t  =  0.63  POP 

R'^  =  0.67 


(1) 


(2) 


where 

t   =  mean  travel  times  in  minutes 

POP  =  urban  area  population. 

Voorhees  et.al.  also  compared  travel  time  variance  with  population.  A 
linear  relationship  involving  travel  time  variance  as  a  function  of 
population  was  determined.   However  the  sample  size  was  small  and 
within  a  narrow  population  range,  and  therefore  the  statistical  strength 
of  the  relationship  was  severely  limited. 

The  main  objective  of  Chan's  research  was  to  synthesize  the 
demand-forecasting  experiences  of  the  past  and  to  come  up  with  an 
estimation  procedure  which  would  allow  the  determination  of  passenger 
travel  on  an  urbanized-area-level. 

Chan  analyzed  the  mean  travel  time  for  the  home  based  work  trips 
with  data  from  86  cities  with  population  ranging  from  55,000  to  4,042,000. 


He  classified  the  urban  areas  by  population  and  structural  form.   For 
those  urban  areas  with  a  population  between  50,000  and  800,000, 
Equations  3  and  4  were  developed  by  regression  analysis  for  a  core- 
concentrated  urban  structure  and  a  multinucleated  urban  structure 
classification. 


Core-concentrated  t  =  6.715  +  0.017  POP 

R^  =  0.48 


(3) 


(4) 


Multinucleated  t  =  -3.748  +  2.134  In  POP 

R^  =  0.83 

where 

t   =  average  travel  time  in  minutes 
POP  =  population  size  in  thousands 
In  =  natural  log  function 

Summary  of  Travel  Time  Statistics  for  the  Study  Areas  in  Indiana 
In  Table  1  are  given  the  location,  base  year  population,  and 
average  travel  time  and  variance  by  trip  purpose  for  each  of  the  study 
areas.   This  Table  clearly  demonstrates  the  wide  range  over  which  the 
areawide  travel  time  statistics  can  vary.   The  minimum  and  maximum 
average  travel  times  for  HBW  trips,  excluding  Indianapolis,  are  8.69 
and  19.71  minutes.   The  minimum  and  maximum  average  travel  times  for 
the  HBO  trips,  excluding  Indianapolis,  are  7.25  and  16.28  minutes. 
Average  travel  times  are  7.37  and  17.00  minutes  for  the  minimum  and 
maximum  values  for  NHB  trips,  excluding  Indianapolis.   Minimum  and 
maximum  travel  time  variances  for  the  HBW,  HBO  and  NHB  trips,  excluding 


en 

(U 

u 

S 

c 

■H 

e 
a 

•H 


0) 
CO 
O 

P. 
M 

3 

a* 

o- 
•1-1 

H 

to 
u 


CO 

•rt 
■U 
CO 
4J 
C/3 

0) 

e 


r-l 

CO 
U 
H 


CO 

CO 
(U 
M 

<: 

13 

3 

4J 

cn 

CO 
•H 

-a 

c 


CO 


PQ 


O 


3 


o 

c 

•H 
>-i 

> 


0) 
60 

0) 
> 


0) 

o 

c 

•H 
U 
CD 

> 


0) 
U 
C 
CO 
•H 
S-i 
CO 
> 


OJD 
CO 
>-. 
0) 
> 


>-• 

C 

CO 

O 

QJ 

•H 

>- 

4-1 

CO 

OJ 

,— 1 

w 

3 

CO 

D. 

P3 

O 

(X, 


4-1 
•H 
O 


O 


c» 


00 


o 


CNl 


o 


o 
in 


CO 


o 
in 


O 

o 


c» 


o 


00 


-J- 

00 


CN 


vD 

00 


CO 


c^ 


CN 

1^ 


O      O 


00 

CTv 

in 


in 
in 


00 

en 


00 

00 


in 


c^     o 
o      i^ 


in 


o 


o 


m     in 


-3- 


in 


in 
t-o 


O 


00     <~n     CXI     in 


o 
o 
o 

00 


O 

o 
o 

CJN 


o 
o 
o 
<:^ 

CO 


o 
o 
o 
o 


o 
o 
o 
o 
o 


o 
o 
o 


O 
O 

o 
o 


O 
O 

o 


en 

3 


O 


•u 
OJ 
>. 

CO 

14-1 

CO 


C 
O 
en 
M 
CD 

-a 
c 
<: 


D 
CO 

a; 
u 

H 


cj 

C 
3 


> 

tn 
c 
CO 

> 


c 
>> 

CO 

3 


O 


13 
C 
QJ 


3 

o 

CO 


O 
O 


00 


ni 

-;; 

-X 

-1^ 

i< 

■K 

nn 

c^ 

^i5 

r~, 

CO 

a^ 

o 

in 

0^ 

vD 

CO 

m 

in 

in 

CM 

ro 

CJ^ 

^ 

CJN 

CJ 

1-1 

. 

• 

• 

• 

• 

■ 

• 

• 

• 

11 

CO 

00 

00 

vD 

r-- 

CTi 

r-~ 

O 

rH 

> 

< 

rH 

1— 1 

r-t 

(TV 


o 

00 


o 

CN 


O 

o 
o 
in 


o 
& 

CO 

c 

CO 

•H 

-o 
c 


CU 
CO 
CO 
Xi 

Q) 

E 
O 

JD 

CO 
3 
O 
•H 
U 
CO 
> 

14-I 
O 

c 

o 


CO 

c 

•H 

B 

o 
o 

-a 

<D 

4-1 

x: 

M 

•H 
01 


T3 

01 
X) 
3 

rH 
CJ 
C 


O 

c 
en 
u 

3 
V^ 

H 

-X 


10 


Indianapolis,  are  9.09  and  42.7  minutes,  9.27  and  40.48  minutes,  and 
6.51  and  50.43  minutes,  respectively. 

Evaluation  of  Travel  Time  Data  for  Small  Urban  Areas 
The  emphasis  of  the  present  research  project  was  on  small  urban 
areas.   For  this  reason,  data  from  Voorhees'  study  and  Chan's  study 
were  combined  into  a  single  data  set  along  with  the  information 
collected  for  this  current  study  for  only  those  urban  areas  with  a 
population  between  38,000  and  400,000  persons.   For  the  sake  of 
uniformity  this  combined  data  set  included  only  home  based  work  trips. 

Figure  1  is  a  plot  of  work  trip  mean  travel  times  versus  popula- 
tion size.   It  can  be  clearly  seen  that  trip  duration  is  positively 
correlated  with  population,  but  there  is  also  a  wide  dispersion  of  the 
mean  travel  times  at  any  population  level.   It  is  interesting  to  note 
that  the  data  points  for  seven  out  of  the  eight  Indiana  urban  areas 
considered  fall  within  a  relatively  narrow  band.   It  is  even  more 
noteworthy  that  one  of  the  Indiana  urban  areas,  Terre  Haute,  is  very 
much  an  outlier.   Linear  regression  analysis  for  the  combined  home 
based  work  data  set  resulted  in  Equation  5.   The  regression  equation 
is  also  plotted  on  Figure  1.   An  equation  was  also  fitted  to  the  seven 
Indiana  urban  areas  that  appear  to  fall  approximately  within  a  narrow 
band  and  is  given  in  Equation  6  for  the  purpose  of  comparison  study. 
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Combined  Data         t=8.4H-0. 01 7P0P     R^=0.23  (5) 

Seven  Indiana  Areas    t=8. 61+0. 018POP     R^=0.33  (6) 

where 

t   =  average  travel  time  in  minutes 

POP  =  population  size  in  thousands 

The  most  important  aspect  of  Figure  1  and  the  regression  equations 
is  not  the  positive  correlation  of  trip  duration  with  population  size, 
but  the  large  spread  of  the  data  points  with  respect  to  the  regression 
line.   Considering  the  above  seven  Indiana  urban  areas  only  and  their 
regression  equation,  average  travel  time  percentage  error  ranged  from 
0.03  percent  to  42.81  percent,  with  the  average  percentage  error  being 
14.20  percent.   This  average  value  is  much  greater  than  FHWA's  suggested 
+  3.0  percent  limit  (13).   Prediction  of  trip  duration  for  home  based 
work  trips  based  on  population  size  alone  would  be  unreliable  due  to 
this  wide  variance  in  the  data. 

Scattergrams  and  regression  analysis  were  also  performed  using 
the  data  collected  in  the  present  study  for  home  based  other  (HBO)  and 
non-home  based  (NHB)  trips.   Positive  correlation  between  mean  travel 
time  and  population  size  was  observed,  but  a  very  wide  dispersion  of 
the  data  points  was  present.   The  large  spread  of  the  data  observed 
severely  limits  the  usefulness  of  any  relationship  developed  for  trip 
duration  prediction  as  was  shown  for  the  above  HBW  trip  data. 

Comparison  of  Travel  Time  Frequency  Distributions 
In  Figures  2  and  3  are  shown  the  travel  time  frequency  distributions 
for  the  home  based  work  trips  for  the  eight  Indiana  study  areas.   The 
division  of  the  eight  areas  into  two  groups  is  for  reasons  of  clarity 
only.   The  Figures  do  illustrate  how  differently  the  travel  time 
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characteristics  can  vary  within  a  given  population  range,  even  though 
home  based  work  (HBW)  trips  are  tlicmselvos   the  least  sensitive  with 
respect  to  travel  time  impedances  as  compared  to  home  based  other  (HBO) 
and  non-home  based  (NHB)  trips. 

Population  range  for  Figure  2  is  from  38,000  persons  (Columbus, 
Indiana)  to  100,000  persons  (Muncle,  Indiana).   In  this  range,  Terre 
Haute  has  travel  time  characteristics  much  different  from  the  other 
urban  areas.   The  population  range  for  Figure  3  Is  176,000  (Evansville, 
Indiana)  to  242,000  persons  (South  Bend,  Indiana).  It  can  be  seen 
that  the  average  travel  time  for  Fort  Wayne,  Indiana  is  much  lower 
than  the  other  two  urban  areas. 

A  close  examination  of  the  travel  time  frequency  distributions 
In  Figures  2  and  3  indicated  that  the  high  average  travel  time  for 
Terre  Haute  and  the  low  average  travel  time  for  Fort  Wayne  might  be 
explained  by  the  magnitude  of  the  zonal  terminal  times  for  these  areas. 

Chapter  Conclusions 
Travel  time  analysis  conducted  in  this  study  revealed  that  the  pre- 
diction of  average  travel  time  based  on  the  population  size  of  the  urban 
area  alone  would  be  of  little  practical  use.   This  conclusion  was 
based  on  the  fact  that  even  though  a  positive  correlation  existed 
between  average  travel  time  and  population  size,  the  variance  of  the 
data  about  the  regression  line  was  substantial  as  measured  by  the  low 
R-squared  values.   This  results  in  a  regression  equation  with  a  confi- 
dence interval  too  large  to  be  of  practical  use,  as  was  evidenced  by  the 
large  average  travel  time  percentage  error  observed  for  the  Indiana  study 
areas . 
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A  comparison  of  home  based  work  travel  time  frequency  distribu- 
tions for  the  Indiana  study  areas  suggested  that  more  careful  selection 
of  zonal  terminal  times  might  enable  better  comparisons  to  be  made 
between  the  various  study  areas  of  Indiana.   However,  travel  time 
frequency  distributions  can  be  compared  only  in  terms  of  general  shape. 
Even  though  two  distributions  may  look  similar,  the  average  travel 
times  can  differ  significantly.   A  difference  of  only  one  minute  can 
easily  result  in  large  discrepancies  in  traffic  assignment.   Using 
travel  time  frequency  distributions  of  one  urban  area  to  compare  with 
the  distributions  generated  from  another  urban  area's  synthetic  trip 
distribution  model  is  therefore  of  little  benefit,  unless  it  is  pre- 
viously known  that  the  average  travel  time  in  the  two  areas  are  very 
close. 

The  above  conclusions  point  out  the  important  fact  that  travel  time 
statistics  of  an  urban  area  are  a  unique  result  of  the  area's  character- 
istics and  are  therefore  difficult  to  compare  with,  or  use  for  pre- 
diction in  another  study  area.   Analyses  conducted  in  Chapters  III  and 
IV  further  support  this  conclusion.   This  uniqueness  between  the  urban 
areas  is  best  exemplified  by  the  descriptions  and  discussions  of  the 
spatial  land  use  distributions  and  network  configurations  of  the 
Lafayette,  Anderson,  and  Muncie   study  areas  given  in  Appendix  A.   At 
first  appearance,  these  three  urban  areas  of  approximately  the  same 
population  size  are  very  similar.   All  three  urban  areas  are  divided 
by  a  river  and  have  a  college  or  university  within  their  boundary. 
However,  a  close  inspection  reveals  that  significant  differences  do 
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exist  in  the  land  use  pattern  and  the  associated  network  configuration, 
and  therefore  the  respective  travel  demand  models  respond  differently 
to  changes  in  the  calibration  parameters. 
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CHAPTER  III 
GRAVITY  MODEL  SENSITIVITY  ANALYSIS 

The  gravity  model  has  been  the  most  widely  used  trip  distribution 
technique  mainly  because  it  is  simple  in  concept  and  the  procedure  is 
well  documented.   Mathematically,  the  gravity  model  is  stated  as 
follows  in  Equation  7: 

P.A.F. .  K.  . 

T   =  -i-J-il-il  (7) 

ij     m 

T.   A.F.  .  K.  . 

j  =  l  J  ^J  ^J 

where: 

T..  =  trips  produced  in  zone  i  and  attracted  to  zone  i, 
ij 

P.   =  trips  produced  by  zone  i, 

A.   =  trips  attracted  by  zone  j,  and 

F. .  =  friction  factor  which  expresses  the  effect  of  spatial 

separation  on  trip  interchanges. 
K. .  =  zone-to-zone  adjustment  factors.   (These  factors  are  ig- 
nored when  synthetic  teclmiques  are  used  because  origin- 
destination  data  is  needed  to  determine  these  values.) 
Friction  factors  can  be  represented  by  a  set  of  values  at  one 
minute  time  increments  or  by  an  equation,  such  as  Equation  8. 

F(t)  =  Ct^'e^''  (8) 

From  equation  8,  two  alternative  forms,  Equations  9  and  10  can  also  be 
formulated: 
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F(t)  =01"  (9) 

and 

F(t)  =  C  e"^*^  (10) 

where : 

F(t)  =  friction  factor  at  time  t, 

C    =  constant  of  proportionality  (in  gravity  model  operation, 
this  constant  is  not  needed  due  to  the  nature  of  the 
model  form) , 

a    =  Alpha,  calibration  parameter, 

3    =  Beta,  calibration  parameter,  and 

e    =  base  of  the  natural  logrithras. 
Equations  8,  9  and  10  are  generally  referred  to  as  the  gamma  function, 
inverse  travel  time  function,  and  the  negative  exponential  function, 
respectively.   Although  in  most  transportation  studies  a  set  of 
friction  factor  values  has  been  used  instead  of  mathematical  equations, 
the  parameters  in  the  equations  enable  transportation  planning 
researchers  to  relate  model  sensitivity  with  respect  to  various 
population  sizes,  urban  structural  forms,  trip  purposes,  and  travel 
time  characteristics  of  various  study  areas. 

The  calibration  of  a  gravity  model  basically  involves  repeated 
adjustments  of  a  set  of  friction  factors  until  the  generated  gravity 
model  travel  time  frequency  distribution  replicates  the  0-D  survey 
travel  time  distribution  within  a  predetermined  error  limit.   Allowable 
error  limit  is  usually  specified  as  a  +  3  percent  error  in  the  gravity 
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model  mean  travel  time  when  compared  to  the  0-D  survey  mean  travel 
time  (13).   Other  gravity  model  accuracy  tests  include  screenline 
checks,  topographical  barrier  checks  and  spiderweb  network  checks. 
The  calibration  procedure,  however,  requires  that  an  0-D  survey  has 
been  performed  in  order  to  conduct  the  accuracy  tests  and  the 
calibration  itself. 

In  order  to  eliminate  the  cost  and  time  involved  in  conducting 
an  0-D  survey,  researchers  and  transportation  planners  have  occasionally 
borrowed  friction  factors  from  other  urban  areas  possessing  similar 
characteristics  and  then  applying  them  to  a  new  study  area.   The 
PLANPAC/BACKPAC  General  Information  (14)  manual  suggests  borrowing 
friction  factors  for  travel  simulation  in  small  urban  areas.   However, 
not  much  work  has  been  done  to  examine  the  sensitivity  of  a  gravity 
model  to  any  possible  errors  that  might  be  introduced  by  the  borrowing 
of  friction  factors  from  a  similar  urban  area  or  by  the  use  of 
standardized  friction  factors  or  other  synthetic  techniques.   Because 
this  is  an  important  aspect  associated  with  synthetic  modelling  of 
transportation  demand,  an  analysis  was  performed  in  the  present  study 
to  evaluate  the  performance  of  the  gravity  model  utilizing  various 
synthetic  techniques. 

The  transportation  study  areas  of  Lafayette  and  Anderson  were 
used  for  the  purpose  of  this  evaluation.   Descriptions  of  these  two 
study  areas  are  given  in  Appendix  A.   The  population  sizes  of  these 
two  study  areas  are  approximately  equal,  but  spatial  activity  distri- 
butions and  network  impedances  are  quite  different  from  each  other. 
These  differences  can  possibly  have  a  profound  effect  on  trip  distribu- 
tion modelling. 
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The  research  discussed  in  this  chapter  consisted  of  four  parts. 
First,  tests  were  conducted  to  measure  the  sensitivity  of  the  gravity 
model  with  respect  to  small  errors  in  the  friction  factors.   Simulta- 
neously, these  tests  revealed  the  sensitivity  of  the  gravity  model 
with  respect  to  the  network  characteristics  and  land  use  activity 
distributions  of  the  individual  study  areas.   Second,  the  variation  of 
friction  factors  with  respect  to  population  size  was  investigated. 
This  investigation  involved  the  examination  of  the  concept  of  similar 
urban  areas  from  which  friction  factors  may  be  borrowed.   Third,  an 
evaluation  was  performed  using  the  standardized  a  parameters,  developed 
by  the  Comsis  Corporation  (3),  in  the  gravity  models  for  the  Lafayette 
and  Anderson  urban  areas.   Fourth,  the  friction  factors  of  Lafayette 
and  Anderson,  Indiana  were  exchanged  or  borrowed  from  each  other  in 
order  to  test  how  borrowed  friction  factors  may  affect  trip  distribu- 
tion modelling  accuracies. 

The  Effect  of  Friction  Factors 

Before  discussing  the  sensitivity  of  gravity  models  with  respect 
to  error  introduced  by  synthetically  determined  or  borrowed  friction 
factors,  it  is  important  to  examine  whether  or  not  friction  factors 
are  necessary  in  trip  distribution  modelling  in  small  urban  areas. 
This  was  accomplished  by  applying  to  the  gravity  model  a  set  of  friction 
factors  of  constant  value  at  all  time  values.   This  results  in  a  trip 
distribution  in  which  relative  travel  times  are  not  a  determinant. 

In  order  to  have  the  friction  factors  equal  a  constant  value,  the 
6  parameter  of  the  negative  exponential  friction  factor  function. 
Equation  10  can  be  set  to  zero.   This  results  in  friction  factors 
being  equal  to  1.0  at  all  time  values.   This  analysis  was  conducted 
with  the  data  from  both  the  study  areas  of  Lafayette  and  Anderson. 
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The  resulting  average  travel  times  as  generated  by  the  gravity 
models  with  constant  friction  factors  are  shown  in  Tables  2  and  3  for 
the  two  urban  areas  considered.   Also  shown  in  the  Tables  are  the  0-D 
survey  mean  travel  times  and  percentage  error  by  trip  purpose.   For 
Lafayette  the  gravity  model  generated  average  travel  time  was  30 
percent  greater  than  the  0-D  survey  value,  while  the  generated  average 
travel  time  for  Anderson  was  19  percent  greater  than  the  corresponding 
0-D  survey  value.   This  relative  difference  in  average  travel  time 
percentage  error  of  11  percent  between  the  two  urban  areas  may  be 
attributable  to  differing  spatial  land  use  distributions  and  network 
characteristics  and  how  the  individual  gravity  model  responds  to 
these  differences. 

The  effect  of  a  constant  friction  factor  was  further  evaluated 
by  comparing  the  travel  time  frequency  distributions  as  determined  by 
the  0-D  survey  and  as  generated  by  the  gravity  model  by  trip  purpose. 
Figures  4  through  9  compare  the  gravity  model  and  0-D  survey  travel 
time  frequency  distributions  for  Lafayette  and  Anderson  by  trip 
purpose.   As  can  be  readily  seen  from  the  Figures,  the  travel  time 
frequency  distributions  of  the  gravity  model  and  0-D  survey  differ 
greatly,  especially  for  the  home  based  other  and  non-home  based  trips. 

The  next  phase  of  the  evaluation  of  the  results  of  a  gravity 
model  with  constant  friction  factors  involved  assignment  of  the  model 
generated  trips  to  the  highway  network  and  subsequent  comparison  of 
the  assigned  trips  to  those  obtained  by  assigning  the  0-D  survey  trips 
to  the  same  highway  network.   The  traffic  assignment  was  based  on  an 
all-or-nothing  assignment  technique  where  the  trips  between  two  zones 
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Table  2.   Results  of  the  Lafayette  Gravity  Model  With  Constant 
Friction  Factor  Values 


Average  Travel  Time 

Trip  Purpose        Survey         Model  Percent  Difference 

HBW  10.57         12.32  +16.56 

HBO  8.56         11.78  +37.62 


0.57 

12.32 

8.56 

11.78 

9.12 

11.42 

9.05 

11.76 

NHB  9.12         11.42  +25.22 

Total  9.05         11.76  +29.94 


Table  3.   Results  of  the  Anderson  Gravity  Model  With  Constant 
Friction  Factor  Values 


Average  Travel  Time 

Trip  Purpose        Survey  Model  Percent  Difference 

HBW             11.44  12.10              +  5.77 

HBO             10.04  12.30              +22.51 

NHB              9.50  11.72              +15.58 

Total            10.14  12.07              +19.03 
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are  assigned  to  the  minimum  time  path  without  regard  to  the  capacities 
of  the  individual  highway  network  links  along  the  path.   An  all-or- 
nothing  assignment  was  used  in  order  to  eliminate  possible  additional 
differences  attributable  to  a  capacity-restrained  assignment  alone  and 
not  to  the  gravity  model. 

The  statistical  test  for  comparing  two  different  traffic  assign- 
ments is  based  upon  grouping  the  individual  links  into  specific  volume 
groups  and  computing  various  statistics  within  each  volume  group.   The 
statistic  most  widely  used  is  the  Percent  Root-Mean-Square-Error  (PRMSE) 
as  defined  in  Equation  11. 


2 
I (test  volume-base  volume) 


PRMSE  =^ ^TT ^^ ^ X  100.0   (11) 

,,     ,         ,          L(base  volume) 
(by  volume  group)  ^t 


where: 

N  =  number  of  links  within  the  base  volume  group. 
An  equation  which  estimates  what  the  expected  PRMSE  should  be  for  a 
given  volume  and  0-D  survey  dwelling  unit  sample  rate  was  developed 
by  A.  Sosslau  and  G.  Brokke  (13) .   A  comparison  can  then  be  made 
between  the  expected  and  observed  values. 

Volume  group  link  comparisons  are  given  in  Tables  4  and  5  for 
Lafayette  and  Anderson,  respectively.   The  0-D  survey  sampling  rate 
for  both  urban  areas  was  10  percent.   It  can  be  seen  that  all  of  the 
observed  PRMSE 's  are  much  larger  than  the  expected  PRMSE  values  for 
the  10  percent  sampling  rate.   This  comparison  is  also  presented 
graphically  in  Figures  10  and  11  for  the  two  urban  areas  considered. 
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Table  4.   All-or-No thing  Traffic  Assignment  Comparison  for  Lafayette 
Gravity  Model  with  Constant  Friction  Factors  Versus  Origin- 
Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

216 

1.48 

34.63 

75.13 

500-  1000 

116 

95.32 

20.25 

51.78 

1000-  2000 

197 

523.15 

14.43 

55.33 

2000-  3000 

140 

930.19 

11.25 

58.22 

3000-  5000 

172 

2023.06 

8.94 

77.39 

5000-10000 

163 

5441.60 

6.54 

98.81 

10000-15000 

31 

10710.32 

5.12 

103.92 

15000-20000 

1 

11281.00 

4.35 

64.90 

Total 

1036 

1759.57 

^— . 

129.83 

Excludes  centroid  links. 


Table  5.   All-or-Nothing  Traffic  Assignment  Comparison  for  Anderson 
Gravity  Model  with  Constant  Friction  Factors  Versus  Origin- 
Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

604 

35.80 

34.63 

84.96 

500-  1000 

211 

240.55 

20.25 

50.63 

1000-  2000 

228 

675.49 

14.43 

56.40 

2000-  3000 

169 

1102.44 

11.25 

53.48 

3000-  5000 

185 

1924.28 

8.94 

54.35 

5000-10000 

108 

3545.44 

6.54 

59.66 

10000-15000 

2 

5440.50 

5.12 

61.60 

Total 

1507 

771.39 

— •" 

86.00 

Excludes  centroid  links. 
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The  equation  for  estimating  the  expected  PRSME  for  a  given  household 
sampling  rate  is  also  shown  in  the  Figures.   It  can  be  seen  from  the 
Figures  that  a  wide  divergence  exists  between  the  expected  and  ob- 
served PRMSE  values  for  both  Lafayette  and  Anderson. 

A  final  test,  a  link  by  link  model/OD  load  ratio  check  was  also 
performed.   This  test  was  performed  on  the  basis  of  the  research  done 
by  Creighton  and  Hamburg  on  the  effect  of  assignment  inaccuracy  on  the 
design  process  in  relation  to  the  number  of  lanes  needed  to  serve  a 
given  hourly  volume  (4).   These  load  ratio  ranges,  as  referenced  in 
the  FHWA  publication.  Traffic  Assignment  (15) , were  used  to  find  the 
number  of  links  with  unacceptable  load  ratios.   The  acceptable  load 
ratio  ranges  used  are  given  in  Table  6. 

Table  6.   Acceptable  Model/OD  Traffic  Assignment  Link  Load  Ratio  Ranges 


Acceptable  Model/OD 
ADT  Volume  Link  Load  Ratio  Range 

1,000  -  5,000  0.55  -  1.45 

5,000  -  10,000  0.60  -  1.40 

10,000  -  15,000  0.65  -  1.35 

15,000  -  20,000  0.70  -  1.30 

20,000  -  30,000  0.75  -  1.25 

Source:  U.S.  Department  of  Transportation,  Federal  Highway  Administra- 
tion, Urban  Planning  Division,  Traffic  Assignment,  Washington, 
D.C. ,  August  1973. 


It  should  be  noted  that  the  link  load  ratios  are  a  result  of  a 
gravity  model  and  an  0-D  survey  all-or-nothing  assignment  and  are  not 
the  loaded  link  volume/ground  count  ratio  normally  used  for  the  test. 
Because  the  purpose  of  the  present  study  is  to  investigate  the 
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sensitivity  of  synthetic  trip  distribution  modelling,  additional  errors 
that  may  be  associated  with  a  capacity-restrained  assignment  and 
errors  in  ground  counts  were  considered  undesirable. 

The  number  of  unacceptable  model/OD  survey  link  load  ratios 
found  for  Lafayette  and  Anderson,  are  respectively  363  and  435  links. 
These  links  represent  35.1  and  28.9  percent  of  all  loaded  links, 
respectively. 

These  results  clearly  point  out  the  importance  of  friction  factors 
in  gravity  model  trip  distribution  simulation  in  small  urban  areas  and 
therefore  also  the  importance  of  travel  time  in  the  trip  making  process. 
It  is  imperative  therefore  that  friction  factor  values  be  calibrated 
as  accurately  as  possible  and  that  it  may  be  necessary  to  exercise 
great  care  when  using  synthetic  methods  in  trip  distribution  modelling. 

Calibration  of  the  Hravlty  Models  for  Lafayet.te  and  Anderson,  Indiana 

The  purpose  of  this  present  study  was  to  evaluate  various 
synthetic  trip  distribution  techniques  and  to  ultimately  recommend  one 
procedure.   Because  the  traditional  gravity  model  calibration  procedure 
is  an  accepted  methodology,  the  trip  distribution  results  from  such 
an  accepted  procedure  were  used  to  serve  as  the  data  base  to  which 
synthetic  methods  were  to  be  compared.   The  objective  was  to  choose  or 
develop  a  teclmiquc  which  can  replicate  the  0-D  survey  trips  almost  as 
accurately  as  the  traditional  method  of  calibration.   It  was  therefore 
necessary  to  (.alibrate  the  Lafayette  and  Anderson  gravity  models  for 
use  in  the  forthcoming  analysis. 
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Because  comparisons  were  to  be  made  between  gravity  models, 
calibration  of  the  gravity  models  for  Lafayette  and  Anderson  using  the 
0-D  data  were  performed  on  the  basis  of  a  negative  exponential  friction 
factor  function.   This  would  allow  comparisons  to  be  made  between 
various  synthetic  techniques  using  the  single  calibration  parameter  6 
instead  of  a  set  of  friction  factor  data  points  which  are  difficult  to 
compare  quantitatively.   Additionally,  the  travel  time  matrices  (skim 
trees)  developed  from  the  minimum  path  trees  and  used  in  the  gravity 
model  calibration  did  not  include  turn  penalties.   Since  the  objective 
of  the  present  research  was  to  simplify  the  trip  distribution  process 
for  small  urban  areas,  it  was  felt  that  the  inclusion  of  turn  penalties 
would  not  appreciably  improve  the  model  results. 

Average  travel  times  and  travel  time  variances,  along  with  the 
3  parameters  obtained  from  the  traditionally  calibrated  gravity  models 
are  given  in  Tables  7  and  8  for  the  study  areas  of  Lafayette  and 
Anderson,  respectively.   Volume  group  link  comparisons,  where  gravity 
model  traffic  assignments  are  compared  to  0-D  survey  traffic  assign- 
ments, were  conducted  for  both  of  the  urban  areas.   Results  of  this 
comparison  are  presented  in  Tables  9  and  10  for  Lafayette  and  Anderson 
respectively.   The  Tables  indicate  that  the  traditionally  calibrated 
gravity  models  resulted  in  very  acceptable  traffic  assignments  for 
both  urban  areas.   Graphical  presentation  of  the  volume  group  link 
comparisons  for  Lafayette  and  Anderson,  are  shown  in  Figures  12  and 
13,  respectively.   A  final  model/0-D  survey  link  load  ratio  test  re- 
vealed that  there  were  no  links  with  an  unacceptable  load  ratio  for 
Lafayette,  while  only  two  links  were  found  to  have  unacceptable  link 
load  ratios  for  Anderson. 
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Table  7.   Traditionnlly  Calibrated  Gravity  Model  Statistics  for 
Lafayette 


Average 

Travel  Time 

Trip 
Purpose 

0-D 
Survey 

Model 

Model  Travel 
Time  Variance 

B 
Parameter 

HBW 

10.60 

10.58 

20.54 

0.1679 

HBO 

8.57 

8.57 

16.11 

0.2856 

NHB 

9.12 

9.12 

12.56 

0.2647 

Total 

9.05 

9.04 

16.18 

— 

Table  8.   Traditionally  Calibrated  Gravity  Model  Statistics  for 
Anderson 


Average 

_T 

ravel  Time 

Trip 
Purpose 

0-D 
Survey 

Model 

Model  Travel 
Time  Variance 

6 
Parameter 

HBW 

11.44 

11.44 

9.09 

0.1327 

HBO 

10.05 

10.04 

8.53 

0.2961 

NHB 

9.50 

9.50 

6.51 

0.3516 

Total 

10.14 

10.14 

8.46 

— 
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Table  9.  All-or-Nothing  Traffic  Assignment  Comparison  for  Lafayette 
with  Traditionally  Calibrated  Gravity  Model  Versus  Origin- 
Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

216 

10.4 

34.63 

49.54 

500-  1000 

116 

32.2 

20.25 

17.13 

1000-  2000 

192 

48.4 

14.43 

13.74 

2000-  3000 

140 

122.0 

11.25 

12.48 

3000-  5000 

172 

176.5 

8.94 

11.40 

5000-10000 

163 

109.0 

6.54 

8.93 

10000-15000 

31 

297.5 

5.12 

4.65 

15000-20000 

1 

-320.0 

4.35 

1.84 

Total 

1036 

86.5 

— 

12.48 

Excludes  centroid  links. 


Table  10.   All-or-Nothing  Traffic  Assignment  Comparison  for  Anderson 
with  Traditionally  Calibrated  Gravity  Model  Versus  Origin- 
Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

604 

7.8 

34.63 

33.43 

500-  1000 

211 

22.6 

20.35 

15.50 

1000-  2000 

228 

54.6 

14.43 

12.11 

2000-  3000 

169 

-16.3 

11.25 

10.31 

3000-  5000 

185 

-36.2 

8.94 

8.99 

5000-10000 

108 

-63.9 

6.54 

7.31 

10000-15000 

2 

-717.0 

5.12 

8.24 

Total 

1507 

3.0 

— 

13.28 

Excludes  centroid  links. 
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The  calibration  of  the  gravity  models  for  Lafayette  and  Anderson 
hence  provided  standards  to  which  synthetic  trip  distribution  models 
could  be  easily  compared.   It  should  be  pointed  out  that  though  both 
of  the  urban  areas  are  of  a  similar  population  size,  there  are  some 
land  use  and  network  characteristics  which  are  distinctively  unique 
to  each  area. 

These  unique  characteristics  resulted  in  dissimilar  travel 
patterns  for  the  two  study  areas.   When  constant  valued  friction 
factors  were  previously  used  on  the  Lafayette  and  Anderson  study  areas 
so  as  to  emphasize  the  importance  of  friction  factors,  the  resulting 
gravity  model  travel  patterns  when  compared  to  their  respective  0-D 
survey  results  demonstrated  significant  differences  between  the  two 
study  areas.   These  differences  between  Lafayette  and  Anderson  do 
provide,  however,  a  wider  statistical  inference  space  over  which  to 
test  the  validity  of  synthetic  trip  distribution  modelling  techniques. 

Previous  Studies  on  the  Use  of  Borrowed  Friction  Factors 
Various  studies  have  been  performed  to  test  the  feasibility  of 
using  borro^^7ed  friction  factors.   The  Federal  Highway  Administration 
in  a  publication  entitled  Urban  Trip  Distribution  Factors  (18) , 
summarized  friction  factor  curve  data  from  a  number  of  transportation 
studies,  specifically  for  the  purpose  of  gathering  the  information  so 
that  friction  factor  curves  could  easily  be  borrowed.   Friction  factor 
curve  data  was  categorized  by  trip  purpose  and  urban  area  population 
size  in  the  FHWA  report.   In  addition,  each  category  of  friction 
factor  curves  was  normalized  into  one  comprehensive  curve. 
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In  order  to  test  the  feasibility  of  borrowinj:',  friction  factor 
curves  from  the  report,  an  experiment  was  conducted  by  the  FHWA  using 
friction  factors  developed  for  Tallahassee,  Florida  and  the  normalized 
friction  factor  values  in  a  test  application  to  the  Sioux  Falls,  South 
Dakota  study  area.   Both  Tallahassee  and  Sioux  Falls  are  in  the  same 
population  range.   On  the  basis  of  a  comparison  of  travel  time  fre- 
quency distributions  and  total  vehicle-miles  traveled,  it  was  con- 
cluded that  the  use  of  borrowed  friction  factors  produced  adequate 
results.   However,  the  available  information  did  not  specify  any 
evaluation  of  traffic  assignment  error.   An  analysis  of  the  traffic 
assignment  would  have  been  a  more  reliable  and  revealing  testing 
methodology. 

In  another  study,  the  Comsis  Corporation  (3)  used  regression 
techniques  to  fit  an  inverse  travel  time  friction  factor  equation  to 
the  FHWA  normalized  friction  factor  curves.   The  a  parameter  is  the 
calibration  parameter  of  the  inverse  travel  time  function.   Much  like 
the  FHWA,  Comsis  Corporation  classified  the  a  parameters  by  trip 
purpose  and  population  size.   A  final  report,  in  which  the  test  case 
results  of  these  standardized  a  parameters  are  reported,  is  not 
published  as  of  this  writing.   However,  these  standardized  a  para- 
meters were  applied  to  the  Lafayette  and  Anderson,  Indiana  study  areas 
as  a  part  of  the  present  study  and  are  reported  in  a  later  section  of 
this  Chapter. 

In  a  comprehensive  report  entitled  Transportation  Planning  For 
Small  Urban  Areas  (5),  Grecco,  et.al.,  reviewed  previous  studies  in 
which  borrowed  friction  factors  had  been  used  in  trip  distribution 
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modelling.   The  conclusion  reached  by  Grecco,  et.al.,  was  that  the 
adoption  of  the  "initial"  friction  factors  for  different  trip  purposes 
from  another  study  area  of  similar  size  is  valid  and  considerably 
simplifies  the  modelling  effort.   The  word  "initial"  in  the  above 
sentence  is  important  for  it  implies  that  adjustments  may  be  necessary 
to  achieve  a  satisfactory  trip  distribution.   Grecco,  et.al. 's  report 
did  not,  however,  mention  the  basis  for  such  a  conclusion. 

Bates  (1),  in  a  synthetic  trip  distribution  modelling  study,  used 
an  inverse  travel  time  friction  factor  equation  to  represent  the 
travel  time  impedance  relationship.   The  inverse  travel  time  function 
is  repeated  in  Equation  12  for  reference. 

F(t)  =  Ct'""  (12) 

where 

F(t)  =  friction  factor  at  time  t, 

C    =  constant, 

a    =  alpha,  calibration  parameter,  and 

e    =  base  of  the  natural  logrlthm. 
Bates  used  the  friction  factor  curves  from  one  city,  Cedartown, 
Georgia,  as  a  base  to  which  adjustments  were  made  according  to  the 
maximum  trip  length  (travel  time)  and  average  trip  length  (travel  time) 
of  the  test  cities.   The  adjustments  consisted  of  altering  the  C 
parameter  of  Equation  12  and  it  was  assumed  that  the  a  parameter  was 
constant  for  all  cities.   Unfortunately,  these  adjustments  of  the  C 
parameter  have  no  effect  on  the  relative  time  values  of  the  friction 
factor  equations.   In  the  gravity  model  equation,  as  repeated  in 
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Equation  13,  the  C  parameter,  which  is  a  constant  value  over  time,  is 
in  both  the  numerator  and  denominator  of  the  equation,  and  therefore 
cancels. 

P.A.F.  .    P.A.C  t~"' 

T   =   ^  -J  ^-1  =  -^J (13) 

ij    I.   A.F..    „  ,  ^  -a  ^^^^ 

■'  1  11    ^  A.C  t 

1 

The  only  parameter  that  can  affect  the  relative  time  values  of  the 
inverse  travel  time  function  is  the  a  parameter. 

Evaluation  of  an  Urban  Area's  Sensitivity  to 
Friction  Factor  Curves 
This  evaluation  was  accomplished  by  observing  the  modelling 
errors  induced  by  small  deviations  from  the  traditionally  calibrated 
3  parameter  in  the  negative  exponential  friction  factor  equation. 
Deviations  of  plus  or  minus  0.05  were  applied  to  the  traditionally 
calibrated  p  parameters  of  each  trip  purpose,  home  based  work,  home 
based  other,  and  non-home  based.   The  resulting  trip  tables  by  purpose 
were  then  added  together  and  loaded  onto  the  network  using  an  all-or- 
nothing  assignment  procedure.   From  these  traffic  assignments  and  the 
0-D  survey  internal-internal  trip  table  network  assignment,  volume 
group  link  comparisons  were  performed.   This  procedure  was  conducted 
on  both  the  Lafayette  and  Anderson  study  areas.   It  was  then  possible 
to  make  relative  comparisons  of  the  magnitude  of  error  introduced  by 
this  evaluation  procedure.   These  results  would  help  to  ascertain 
whether  different  synthetic  techniques  may  be  necessary  for  different 
urban  areas  or  if  a  single  synthetic  technique  can  be  found  to 
adequately  model  the  travel  patterns  of  the  uniquely  defined  urban  areas. 
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Lafayette's  gravity  model  average  travel  time  statistics 
demonstrating  S  parameter  deviation  sensitivity  are  given  in  Table 
11.   Percentage  change  in  average  travel  times  were  approximately  of 
the  same  magnitude  for  both  the  plus  and  minus  0.05  deviations  from 
the  traditionally  calibrated  3  parameter.   The  percentage  differences 
between  the  0-D  survey  average  travel  time  and  the  6  parameter  +0'05 
deviations  from  the  traditionally  calibrated  B  parameter  are  greater 
than  the  +3.0  percent  level  suggested  by  the  FHWA  (13). 


Table  11.  Lafayette  Total  Trip  Average  Travel  Time  for  6  Parameter 
Deviations  from  the  Traditionally  Calibrated  6  Value  for 
Lafayette 


3  Parameter  Deviation 
-0.05         0.00       +0.05       +0.10 


Total  Trip  Average 

Travel  Time  9.48         9.05        8.63        8.31 

Percent  Error  +4.75  —        -4.64       -8.18 


A  greater  indication  of  error  as  a  result  of  the  3  parameter 
deviations  from  the  traditionally  calibrated  value  can  be  found  in 
the  volume  group  link  comparisons  as  shown  in  Table  12.   Total  PRJ-ISE 
for  the  minus  0.05  3  deviation  was  almost  twice  that  of  the  traditionally 
calibrated  gravity  model,  while  the  total  PRMSE  of  the  plus  0.05  3 
parameter  deviation  was  only  one  third  greater.   As  a  further  test,  a 
plus  0.10  3  parameter  deviation  was  applied  to  the  Lafayette  data  set. 
The  resulting  total  PRMSE,  also  shown  in  Table  12,  was  a  little  over 
twice  the  normally  calibrated  gravity  model  PRMSE  and  is  very  close 
to  the  minus  0.05  3  parameter  deviation  PRMSE. 
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Table  12.   All-or-Nothing  Traffic  Assignment  Comparison  for  Lafayette: 
3  Parameter  Deviation  Sensitivity  Analysis  * 


Volume 

Expected 

3  -0.05 

Traditional 

3  +0.05 

3  +0.10 

Group 

PRMSE 

Deviate 

Model 

Deviate 

Deviate 

X  1000 

10%  Sample 

PRMSE 

PRMSE 

PRMSE 

PRMSE 

0-  0.5 

34.63 

46.65 

49.54 

53.99 

59.17 

0.5-  1.0 

20.25 

21.93 

17.13 

16.91 

19.79 

1.0-  2.0 

14.43 

18.01 

13.74 

13.60 

16.22 

2.0-  3.0 

11.25 

16.62 

12.48 

12.40 

14.93 

3.0-  5.0 

8.94 

18,82 

11.40 

11.55 

15.97 

5.0-10.0 

6.54 

17.90 

8.93 

12.79 

20.54 

10.0-15.0 

5.12 

17.92 

4.65 

9.73 

18.85 

15.0-20.0 

4.35 

8.11 

1.84 

10.46 

17.53 

Total 



24.45 

12.48 

16.58 

26.45 

'"'Excludes  centroid  links. 

Continuing  this  analysis,  a  link  by  link  comparison  was 
conducted  in  order  to  find  unacceptable  model/O-D  link  load 
ratios  on  the  basis  of  the  0-D  survey  traffic  assignment  and  the 
assignment  of  trip  interchanges  resulting  from  the  different  levels 
of  the  3  parameters.   The  results  of  this  link  load  ratio  analysis 
are  given  in  Table  13.   The  Table  shows  that  the  minus  0.05  3  para- 
meter deviate  assignment  has  more  unacceptable  link  load  ratios  than 
either  of  the  plus  0.05  or  0.10  3  deviates. 


Table  13.   Number  of  Unacceptable  Model/O-D  link  Load  Ratios  for 
Lafayette  3  Deviate  Model  Traffic  Assignments 


3  -0.05 

Deviate 

Model 


Traditionally 

Calibrated 

Gravity 

Model 


3  +0.05 

Deviate 

Model 


3  +0.10 

Deviate 

Model 


Number  of  Un- 
acceptable 
Load  Ratios 


15 
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This  same  3  parameter  deviation  sensitivity  analysis  was  also 
conducted  on  the  Anderson  study  area  data.   Anderson  average  travel 
time  statistics,  demonstrating  the  sensitivity  of  the  gravity  model  to 
3  parameter  deviations  from  the  traditionally  calibrated  3  parameter, 
are  given  in  Table  14.   The  percentage  difference  in  total  trip 
average  travel  times,  when  compared  to  those  obtained  from  the  0-D 
survey,  were  approximately  of  the  same  magnitude  for  both  the  plus  and 
minus  0.05  deviations  in  the  3  parameter.   These  errors  were  approxi- 
mately within  the  allowable  range  as  suggested  by  the  FHWA  (13) . 


Table  14.   Total  Trip  Average  Travel  Time  for  3  Parameter  Deviations 
from  the  Traditionally  Calibrated  3  Parameter  for  Anderson 


3  Parameter  Deviation 
-0.05  0.00  +0.05 

Total  Trip  Average 

Travel  Time  10.47  10.14  9.85 

Percent  Error  +  3.25  —  -2.86 


Volume  group  link  comparison   PRMSE  results   are  shovm  in 
Table  15.   The  total  PRMSE  for  the  plus  and  minus  3  parameter  deviations 
are  very  close,  being  only  about  one  third  more  than  the  traditionally 
calibrated  gravity  model  value. 

A  link,  by  link  comparison  was  performed  in  order  to  find  un- 
acceptable gravity  model  0-D  survey  link  load  ratios  for  the  0-D  survey 
assignment  and  the  assignments  resulting  from  various  3  deviates. 
The  results  of  this  analysis  are  given  in  Table  16.   As  indicated  in 
the  Table,  the  plus  3  parameter  deviate  had  fewer  links  with  unacceptable 
load  ratios  than  the  minus  0.05  3  parameter  deviate. 
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Table  15.   All-or-Nothing  Traffic  Assignment  Comparison  for  Anderson: 
3  Parameter  Sensitivity  Analysis 


Traditionally 

Volume 

Expected 

B  -0.05 

Calibrated 

3  +0.05 

Group 

PRflSE 

Deviate 

Model 

Deviate 

X  1000 

10%  Sample 

PRMSE 

PRMSE 

PRMSE 

0-  0.5 

3A.63 

37.48 

33.43 

34.26 

0.5-  1.0 

20.25 

18.28 

15.50 

16.34 

1.0-  2.0 

14.43 

16.46 

12.11 

15.86 

2.0-  3.0 

11.25 

12.72 

10.31 

13.19 

3.0-  5.0 

8.94 

12.24 

8.99 

11.97 

5.0-10.0 

6.54 

11.13 

7.31 

11.25 

10.0-15.0 

5.12 

2.94 

8.24 

17.62 

Total 

— 

18.28 

13.28 

18.18 

Excludes  centroid  links. 


Table  16.   Number  of  Unacceptable  Model/OD  Link  Load  Ratios  for 
Anderson  S  Deviate  Model  Traffic  Assignment 


3 

-0.05 

Deviate 

Model 


Traditionally 
Calibrated 

3 
+0.05 

Gravity 

Model 

Deviate 

Model 

Number  of  Unacceptable 
Link  Load  Ratios 


The  relatively  low  number  of  unacceptable  links  for  the  Lafayette  and 
Anderson  3  parameter  deviates  are  probably  trivial  because  they  represent 
a  very  small  percentage  of  all  links  and  are  only  one  of  the  many  sources 
of  error  in  travel  demand  modeling.   It  should  be  remembered  however  that 
the  end  result  of  the  entire  research  effort  is  the  implementation  of  a 
completely  synthetic  travel  demand  modelling  system  that  is  capable  of 
modelling  urban  areas  whose  internal  characteristics  may  vary  greatly. 
Therefore,  error  must  be  kept  to  a  minimum  within  each  step  of  the  process 
so  as  not  to  be  carried  into  successive  steps  if  the  resources  spent  to 
obtain  such  accuracy  are  minimal. 
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The  above  analysis  on  the  Lafayette  and  Anderson  study  areas 
suggests  that  the  sensitivity  of  a  gravity  model  with  respect  to 
friction  factor  errors  may  be  dependent  upon  the  particular  urban  area 
that  is  being  modoU'd .   The  errors  introduced  by  the  (■'  parameter 
deviations  from  the  calibrated  value  were  much  greater  for  Lafayette 
than  for  Anderson.   This  varying  sensitivity  may  cause  difficulties  in 
the  selection  of  a  synthetic  trip  distribution  modelling  technique 
applicable  to  all  small  urban  areas. 

Variability  of  Friction  Factor  Data  Within  Population  Size 
The  main  criterion  for  the  selection  of  a  similar  urban  area  from 
which  to  borrow  friction  factors  has  been  the  population  size.   In 
order  to  test  the  validity  of  this  selection  criterion,  a  negative 
exponential  function  was  fitted  to  the  friction  factor  curves  of  various 
urban  areas  collected  for  the  present  study.   This  procedure  provided 
the  S  calibration  parameter  on  which  the  criterion  v\ras  to  be  tested. 
Appendix  B  discusses  the  regression  technique  used  in  determining  the 
3  parameters. 

In  Figure  14,  the  computed  P.  parameters  are  plotted  against 
population  size  for  home  based  work  trips.   The  Figure  shows  that,  for 
this  data  set,  population  size  is  a  poor  criterion  by  which  to  select 
a  similar  urban  area,  primarily  due  to  the  wide  dispersion  of  the  data 

points  at  the  smaller  population  sizes.   The  intercept,  slope,  coeffi- 

2 
cient  of  multiple  determination  (R  ),  and  the  significance  level,  as 

determined  by  linear  regression  analysis  are  also  shown  on  the  Figure. 

Similar  results  were  obtained  for  the  home  based  other  and  non-home 

based  trips  as  shown  in  Figures  15  and  16. 
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Of  partinilar  interest  is  Figure  14  where  the  slope  of  the 
regression  line  is  positive  for  this  data  set.   The  usual  assumption 
is  that  the  p.  parameter  (slope)  decreases  with  increasing  population 
size  (7).   The  positive  slope  can  probably  be  attributed  to  a  small 
sample  size  and/or  the  limited  population  range  to  which  this  line  was 
fitted.   The  assumption  that  the  slope  should  be  negative  implies  that 
the  variability  of  the  B  parameter  with  respect  to  population  size  is 
even  greater  than  what  Figure  14  implies.   Consequently,  it  would 
appear  that  the  borrowing  or  the  use  of  standardized  friction  factors 
or  parameters  may  possibly  result  in  unacceptable  trip  distribution 
models. 

Evaluation  of  Standardized  Friction  Factor  Parameters 
The  ability  to  use  standardized  friction  factors  or  parameters  of 
friction  factor  functions  could  greatly  reduce  the  time  and  resources 
expended  in  gravity  model  calibration.   A  test  was  conducted  by  running 
the  gravity  models  for  Lafayette  and  Anderson  using  the  standardized 
a  parameters  of  the  inverse  travel  time  function.   These  a  parameters 
were  developed  by  the  Comsis  Corporation  using  the  data  contained 
within  FHWA's  publication  entitled  Urban  Trip  Distribution  Factors  (18), 

Lafayette's  gravity  model  travel  time  statistics  which  resulted 
from  the  application  of  standardized  a  parameters,  are  presented  in 
Table  17.   Percentage  differences  of  the  gravity  model  average  travel 
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t Lmcs  from  those  obtaincnl  from  the  O-D  survey  values  are  minor,  the 

largest  value  being  only  a  -2.08  percent  for  the  home  based  work  trips. 

The  standardized  a  parameters  were  also  very  close  to  the  actual 

calibrated  parameters  as  shown  in  Table  17.   The  PRMSE  statistics  as 

determined  from  a  volume  group  link  comparison  are  reported  in  Table  li 

and  also  graphically  illustrated  in  Figure  17.   It  is  evident  from  the 

Table  and  Figure  that  the  PRMSE  values  at  each  volume  level  are  very 

acceptable.   A  link  by  link  check  found  12  out  of  1036  links  with  an 

acceptable  link  load  ratio  as  compared  to  no  links  for  the  traditional 

gravity  model  calibration. 

Table  18.   All-or-Nothing  Traffic  Assignment  Comparison  for  Lafayette: 
Standardized  a  Parameters  Versus  Origin-Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

216 

7.37 

34.63 

47.57 

500-  1000 

116 

44.41 

■  20.25 

23.90 

1000-  2000 

197 

54.38 

14.43 

14.71 

2000-  3000 

140 

64.45 

11.25 

10.97 

3000-  5000 

172 

48.12 

8.94 

10.87 

5000-10000 

163 

-285.23 

6.54 

9.37 

10000-15000 

31 

-349.19 

5.12 

4.27 

15000-20000 

1 

-1499.00 

4.35 

8.62 

Total 

1036 

-23.22 

— 

12.68 

Excludes  centroid  links. 


Gravity  model  travel  time  statistics  for  the  Anderson  study  area, 
as  generated  by  the  inclusion  of  the  standardized  «  parameters,  are 
presented  in  Table  19  along  with  the  O-D  survey  values  for  comparison. 
Percentage  differences  in  model  average  travel  times  from  O-D  survey  values 
bv  trip  purpose  are  relatively  low  in  comparison  to  the  large  absolute 
differences  between  the  standardized  't  parameters  and  the  traditionally 
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calibrated  values.   For  NHB  trips,  the  predicted  average  travel  time 
was  5.74  percent  greater  than  the  0-D  survey  average  travel  time. 
This  error  is  almost  twice  the  FHWA  suggested  +3.0  percent  error  (13) 
The  HBW  trip  error,  however,  was  of  opposite  sign  and  therefore 
helped  to  counteract  the  NHB  trip  error.   Traffic  assignment  volume 
group  link  comparisons  revealed  acceptable  PRMSE  values  at  all  volume 
ranges  as  shown  in  Table  20  and  Figure  18.   Only  three  of  the  1507 
links  were  found  with  unacceptable  model/0-D  survey  link  load  ratios, 
only  one  more  than  the  normally  calibrated  gravity  model. 

Table  20.   All-or-Nothing  Traffic  Assignment  Comparison  for  Anderson: 
Standardized  a   Parameter  Versus  Origin- Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

604 

14.46 

34.63 

37.83 

500-  1000 

211 

64.42 

20.25 

18.85 

1000-  2000 

228 

133.44 

14.43 

15.58 

2000-  3000 

169 

113.65 

11.25 

11.64 

3000-  5000 

185 

176.89 

8.94 

10.65 

5000-10000 

108 

354.29 

6.54 

9.45 

10000-15000 

2 

14.00 

5.12 

1.15 

Total 

1507 

94.87 

— 

16.14 

Excludes  centroid  links. 


The  preceding  analysis  of  using  standardized  a   parameters  on  the 
two  test  cases  help  to  possibly  explain  v>;hy  transportation  planners 
and  researchers  have  usually  concluded  that  in  general,  borrowed  or 
standardized  friction  factors  often  produce  acceptable  synthetic  trip 
distribution  model  results.   Total  trip  as  well  as  the  individual  trip 
purpose  average  travel  time  statistic  and  the  volume  group  link 
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comparison  for  Lafayette  do  indicate  acceptable  model  results.   These 
results  are  duo  to  the  fact  that  the  standardized  and  the  actually 
calibrated  a   parameters  are  very  close  in  magnitude.   The  Lafayette 
results  therefore  suggest  that  the  use  of  standardized  a  parameters 
is  indeed  a  valid  synthetic  procedure  for  trip  distribution  modelling. 
Similarly,  total  trip  average  travel  time  and  volume  group  link 
comparisons  for  Anderson  also  indicate  that  an  acceptable  trip  distri- 
bution model  result  had  occurred  and  therefore  additionally  endorsing 
the  use  of  standardized  a  parameters. 

However,  just  because  the  aggregated  results  appear  acceptable,  it 
cannot  be  necessarily  inferred  that  the  intermediate  steps  to  achieve 
the  final  results  are  accurate  and  completely  acceptable.   Closer 
inspection  of  the  intermediate  Anderson  model  results  reveal  that  NHB 
trip  average  travel  time  has  been  overestimated  by  5.74  percent  when 
compared  to  the  0-D  survey  data.   One  reason  why  this  overestimate 
was  not  detrimental  to  tlie  final  results  was  the  HBW  trip  average 
travel  time  was  underestimated  by  2.29  percent  and  therefore  helped 
to  offset  the  NHB  trip  error.   If  by  chance,  as  it  may  possibly  occur 
for  another  study  area,  the  HBW  trip  average  travel  time  was  also 
overestimated,  then  the  total  error  may  have  been  much  greater.  Another 
possible  reason  for  Anderson's  total  model  results  to  be  acceptable,  is 
the  fact  that  its  gravity  model  sensitivity  to  friction  factor  error 
is  small  as  was  demonstrated  in  the  previous  section  of  this  Chapter. 
Table  19  showed  that  there  is  indeed  a  large  difference  between  the 
standardized  a  parameters  and  the  a   parameters  that  were  determined 
from  traditional  calibration  methods  for  the  HBW  and  NHB  work  trips. 


61 


Evaluation  of  Borrowed  Friction  Factors 
As  a  final  example  of  the  sensitivity  of  gravity  models  with 
respect  to  borrowed  friction  factors,  the  Lafayette  and  Anderson 
gravity  models  were  run  using  each  other's  calibrated  3  parameters 
of  the  negative  exponential  friction  factor  equation.   Both  of  the 
urban  areas  are  approximately  of  the  same  population  size.   The 
generated  trip  tables  were  then  assigned  to  their  respective  networks 
and  results  analyzed. 

Average  travel  time  for  Lafayette's  total  trips  resulting  from 
the  gravity  model  with  Anderson's  borrowed  friction  factor  6  parameter 
was  only  2.1  percent  less  than  that  obtained  from  the  0-D  survey,  as 
given  in  Table  21.   However,  both  the  home  based  work  trips  and  the 
non-home  based  trips  average  travel  times  were  greater  than  the 
allowable  +3.0  percent  error.   These  values,  however,  are  of  opposite 
sign  and  therefore  have  a  tendency  to  cancel  each  other's  error. 
Volume  group  link  comparison  PRMSE  values  reported  in  Table  22  and 
illustrated  in  Figure  19,  are  actually  slightly  better  than  the 
supposedly  more  accurate  traditionally  calibrated  model  values.   A 
check  of  the  model/0-D  survey  link  load  ratios  revealed  no  links  with 
an  unacceptable  load  ratio. 

The  good  performance  of  the  traffic  assignment  is  probably  due 
to  the  offsetting  errors  as  discussed  above.   If  the  home  based  work 
trip  average  travel  time  had  actually  been  underpredicted,  the  traffic 
assignment  would  have  probably  shown  increased  error. 
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Table  21.   Gravity  Model  Statistics  for  Lafayette  with  Anderson  S 
Parameter 


Trip  Purpose 


HBW 
HBO 
NHB 
Total 


Average  Travel  Time 
Survey        Model 


10.60 

10.94 

8.57 

8.47 

9.12 

8.50 

9.05 

8.86 

Percent  Difference 


+  3.21 
-1.17 
-6.80 
-2.10 


Table  22.   All-or-Nothing  Traffic  Assignment  Comparison  for  Lafayette 
with  Anderson  L^.  Parameter  Versus  Origin-Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

216 

6.95 

34.63 

49.20 

500-  1000 

116 

13.53 

20.25 

16.11 

1000-  2000 

197 

10.59 

14.43 

13.09 

2000-  3000 

140 

57.44 

11.25 

11.93 

3000-  5000 

172 

49.96 

8.94 

10.34 

5000-10000 

163 

-210.02 

6.54 

9.30 

10000-15000 

31 

-248.45 

5.12 

3.95 

15000-20000 

1 

-1160.00 

4.35 

6.67 

Total 

1036 

— 

12.31 

Excludes  centroid  links. 
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A  gravity  model  run  for  Anderson,  using  the  (-'-  parameters 
developed  for  Lafayette,  resulted  in  an  average  travel  time  for  all 
trips  of  only  1.8  percent  greater  than  the  corresponding  mean  travel 
time  obtained  from  the  0-D  survey.   Table  23  shows  that  the  non-home 
based  trip  average  travel  time  was  overestimated  by  5.68  percent. 
Home  based  work  trip  average  travel  time  did  however,  help  to 
compensate  for  the  error  of  the  non-home  based  trips.   Volume  group 
link  comparisons,  as  indicated  in  Table  24  and  Figure  20  were 
acceptable.   Only  three  links  out  of  1507  were  found  unacceptable, 
one  more  than  the  traditionally  calibrated  model. 

For  the  two  study  areas  of  Lafayette  and  Anderson,  the  borrowing 
of  each  other's  friction  factor  parameters  indicated  good  results  for 
several  reasons.   First,  Lafayette's  large  percentage  error  in 
prediction  of  average  travel  time  for  two  of  the  three  trip  purposes 
were  of  opposite  sign  and  thus  the  aggregated  total  error  was  within 
the  acceptable  level.   The  same  situation  of  offsetting  errors  took 
place  in  Anderson,  but  of  a  smaller  magnitude.   Secondly, 
Anderson's  gravity  model  is  less  sensitive  to  friction  factor  errors 
as  compared  to  Lafayette. 

Chapter  Conclusions 
Sensitivity  analysis  of  the  gravity  models  for  Lafayette  and 
Anderson  resulted  in  three  important  observations  in  regards  to 
synthetic  trip  distribution  modelling. 

1.   Acceptable  model  results  will  occur  when  the  borrowed,  or 

standardized  friction  factors  are  very  similar  to  those  that 
would  have  been  normally  calibrated  from  0-D  survey  data. 
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Table  23.   Gravity  Model  Statistics  for  Anderson  with  Lafayette  3 
Parameter 


Trip  Purpose 


Average  Travel  Time 
Survey         Model 


Percent  Difference 


HBW 
HBO 
NHB 
Total 


11.44 

11.25 

10.05 

10.11 

9.50 

10.04 

10.14 

10.32 

-1.66 

+0.60 
+5.68 
+1.78 


Table  24.   All-or-Nothing  Traffic  Assignment  Comparison  for  Anderson 

with  Lafayette  B  Parameter  Versus  Origin-Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRIISE 

PRMSE 

0-   500 

604 

13.71 

34.63 

36.01 

500-  1000 

211 

50.87 

20.25 

17.20 

1000-  2000 

228 

114.05 

14.43 

14.37 

2000-  3000 

169 

84.30 

11.25 

10.70 

3000-  5000 

185 

125.01 

8.94 

9.67 

5000-10000 

108 

217.92 

6.54 

8.17 

10000-15000 

2 

-269.50 

5.12 

3.35 

Total 

1507 

69.93 

— 

14.48 

Excludes  centroid  links. 
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2.  Acceptable  model  results  may  occur  when  synthetic  methods 
are  used  because  the  urban  area's  gravity  model  is  not 
sensitive  to  friction  factor  errors  (as  was  the  case  for 
Anderson) . 

3.  Acceptable  model  results  may  occur  when  one  individual  trip 
purpose  error  is  offset  by  an  opposite  error  in  one  or  more 
of  the  other  trip  purposes. 

Clearly,  the  first  situation  is  the  most  preferred  with  the 
second  also  being  acceptable.   The  third  situation,  however,  is  un- 
acceptable because  chance  plays  a  larger  role  in  the  model's  outcome 
than  the  procedure  by  which  the  friction  factors  were  obtained. 

In  addition,  analysis  of  data  from  a  number  of  small  urban  area 
transportation  studies  resulted  in  one  other  important  conclusion  in 
regards  to  the  borrowing  of  friction  factors.   Population  size  is  an 
inadequate  criterion  by  which  to  select  a  similar  urban  area  for 
borrowing  friction  factors  or  their  parameters. 

The  above  observations  and  conclusions  suggest  that  in  many 
instances,  the  borrowing  of  friction  factors  or  the  use  of  standardized 
values  in  trip  distribution  gravity  models  may  result  in  total  trip 
models  of  undesirable  accuracy.   Therefore,  a  technique  which  produces 
consistently  more  reliable  results  is  still  desired,  for  again,  it 
must  be  remembered  that  synthetic  trip  distribution  modelling  is  only 
one  step  within  a  sequence  of  models. 
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CHAPTER  IV 
DEVELOPMENT  AND  TESTING  OF  PROCEDURES  FOR  SYNTHETIC 
SELF- CALIBRATION  OF  GRAVITY  MODELS 

In  the  previous  Chapter,  the  two  major  methods  of  gravity  model 
calibration  -  borrowed  friction  factor  curves  or  their  parameters   and 
the  use  of  standardized  friction  factor  values  were  evaluated.   The 
results  indicated  that  these  methods  can  produce,  in  some  cases,  in- 
consistent and  unreliable  results.   This  discrepancy  in  the  results 
can  be  attributed  primarily  to  the  variability  of  various  urban 
characteristics  among  the  many  small  urban  areas.   Thus,  there  exists 
a  definite  need  to  have  a  self-calibrating  synthetic  gravity  model  which 
would  be  reliable  in  reflecting  the  varying  characteristics  of  different 
urban  areas.   Such  a  model  would  significantly  reduce  the  time  and 
effort  presently  expended  in  the  development  of  a  trip  distribution 
model  for  an  urban  area  without  introducing  any  significant  error  into 

the  modelling  process. 

Literature  Review 

Most  of  the  research  on  self-calibrating  gravity  models  has 

centered  on  the  approximation  of  the  3  parameter  of  the  negative 

exponential  friction  factor  equation.   Researchers  such  as  Tomlin  (12), 

Wardrop  (20) ,  and  Zaryouni  (22)  suggested  the  approximation  of  6  as 

the  inverse  of  the  areawide  mean  travel  time.  Equation  14,  without 

regard  to  any  particular  trip  purpose. 
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3  :^  1/C  (1^) 

where 

C  equals  areawide  average  travel  time. 
An  equation  of  this  type  can  easily  be  incorporated  into  an  iterative 
procedure  for  calibrating  a  gravity  model.   The  testing  of  such  a 
procedure,  however,  has  been  minimal  with  both  favorable  and  unfavorable 
results. 

Another  possible  approach  to  synthetic  self-calibration  of  the 
gravity  model  may  be  at  a  stratified  level.   Voorhees  (19),  in  an 
effort  to  measure  the  spatial  arrangement  of  activities,  introduced 
the  concept  of  the  opportunity  distribution;  that  is,  the  gravity 
model  travel  time  distribution  resulting  from  the  setting  of  the 
friction  factors  equal  to  a  constant  value  of  1-0.   Voorhees  used  the 
1948  Washington,  D.  C.  transportation  data  to  reveal  how  differently 
opportunity  travel  time  frequency  distributions  can  vary  between  zones. 
It  was  concluded  that   "in  the  case  of  the  gravity  model,  separate 
travel  time  factors  should  be  used  if  there  is  a  wide  variation  in 
opportunity  distribution  between  zones".   Consequently,  Voorhees 
suggested  the  development  of  a  new  distribution  model  that  recognizes 
both  travel  time  and  spatial  arrangement  of  opportunities. 

Further  support  of  this  possible  approach  was  provided  in  the 
article  by  Fisk  and  Brown  (5)  which  states:  "The  parameter  8  [of  the 
negative  exponential  friction  factor  equation]  was  found  to  be  sensi- 
tive to  the  distribution  of  destinations  about  the  origin  zone...". 
Wilson  (21)  also  suggested  that  the  3  parameters  should  be  origin 
specific. 
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The  literature  review  thus  provides  two  possible  approaches  to 
synthetic  self-calibration  of  the  gravity  model;  1)  areawide  parameter 
estimation  using  areawide  travel  time  statistics,  and  2)  origin  zone 
specific  parameter  estimation  using  origin  zone  travel  time  statistics. 

Areawide  Approach  to  Synthetic  Self-Calibration  of  a  Gravity  Model 
The  friction  factor  equation  of  a  gravity  model  represents  the 
trip  makers'  destination  choice  process  with  respect  to  the  spatial 
distribution  of  the  land  use  activities  as  represented  by  the 
destination  choice  variable  of  interzonal  travel  time.   The  aggregate 
effect  of  all  the  trip  makers'  destination  choices  result  in  a  total 
travel  pattern  that  can  be  approximately  described  in  terms  of  the 
areawide  travel  time  statistics  of  mean  and  variance.   Average  travel 
time  and  travel  time  variance  can  thus  be  thought  of  as  areawide 
descriptors  of  the  entire  distribution  of  land  use  activities  and  the 
type  and  configuration  of  the  network.   It  can  be  expected,  therefore, 
that  it  should  be  possible  to  relate  mathematically  the  parameter  of 
a  friction  factor  equation  to  the  spatial  land  use  activity  distribu- 
tion and  network  impedances  through  the  use  of  the  areawide  descriptors 
of  mean  travel  time  and  travel  time  variance,  and  possibly  to  other 
variables  of  interest.   These  relationships  are  represented  by  the 
schematic  diagram  in  Figure  21. 

Data  used  in  the  study  of  the  areawide  approach  consisted  of  the 
B  parameters  of  the  negative  exponential  friction  factor  function,  area- 
wide  average  travel  times  and  travel  time  variances,  study  area  population, 
and  study  area  density  as  computed  or  reported  for  eighteen  small  urban 
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area  transportation  studies.   Regression  analysis  was  used  to  search 
for  mathematical  relationships  between  the  R  parameters  and  the 
areawide  travel  time  statistics  and  any  other  variables  of  possible 
interest  in  the  first  approach. 

Extensive  linear  and  nonlinear  regression  analysis  resulted  in 
mathematical  relationships  of  very  simple  form  for  each  trip  purpose. 
Equations  15,  16,  and  17  give  the  relationship  between  the  3  parameters 
of  the  negative  exponential  friction  factor  function  and  average 
travel  time  for  the  home  based  work  (HBW) ,  home  based  other  (HBO) ,  and 
non-home  based  (NHB)  trips  respectively.   The  coefficient  of  multiple 
determination,  R   is  also  given. 


HBW:   6  =  0.0594  +  0.7858/C  (15) 

R^  =  .278 


HBO:   3  =  0.01A5  +  2.0794/C  (16) 

R^  =  .479 


NHB:   6  =  -0.0222  +  2.2976/C  (17) 

2 
where  R  =  . 502 

C  =  areawide  average  travel  time. 


The  forms  of  Equations  15,  16,  and  17  are  similar  to  that  of 
Equation  14.   However,  the  HBW  equation  (Equation  15)  contains  a 
significant  constant  term  while  the  HBO  and  NHB  equations  have 
insignificant  constant  terms.   Furthermore,  the  slopes  (coefficients) 
of  the  independent  variable  vary  significantly  with  trip  purpose. 
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Equation  15  for  HBW  trips  is  the  only  equation  whose  slope  (coefficient) 
is  close  to  that  of  hypothesized  form  Equation  14.   The  slopes 
(coefficients)  do,  however,  emphasize  the  varying  degree  of  importance 
placed  on  travel  time  for  the  three  trip  purposes. 

The  statistical  fit  of  the  Equations  are  less  than  what  would 
normally  be  considered  acceptable  for  estimating  purposes.   Neverthe- 
less, because  other  researchers  have  suggested  that  relationships  of 
this  form  may  produce  acceptable  model  results,  these  Equations  were 
incorporated  into  an  iterative  algorithm  for  calibrating  a  gravity 
model  so  as  to  evaluate  this  technique.   A  simplified  flow  chart  of 
this  iterative  procedure  is  given  in  Figure  22.   Because  of  data 
availability  and  the  previously  noted  differences  in  each  urban  area's 
model  sensitivity  to  friction  factor  errors,  the  Lafayette  and 
Anderson  study  areas  were  selected  for  validation  testing. 

Tables  25  and  26  contain  the  average  travel  time  statistics  by 
trip  purpose,  as  generated  by  the  test  gravity  model  for  the  urban 
areas  of  Lafayette  and  Anderson,  respectively.   Also  sho;\m  in  Tables 
25  and  26  are  the  comparisons  of  these  average  travel  time  values  with 
those  obtained  from  the  origin-destination  surveys.   Total  trip 
average  travel  times  were  overestimated  by  the  areawide  test  gravity 
model  for  both  Lafayette  and  Anderson  by  3.76  percent  and  6.35  percent, 
respectively.   In  addition,  the  3  parameters  as  computed  by  the  test 
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INITIAL  ESTIMATE  OF  FRICTION 
FACTOR  3  PARAl^lETER 


DISTRIBUTE  TRIPS 


COMPUTE  THE  RESULTING  AREAWIDE 
AVERAGE  TRAVEL  TIME 


COMPUTE  NEW  ITERATED  g  PARAMETER 

AS  A  FUNCTION  OF  AREAWIDE 

AVERAGE  TRAVEL  TIME 


CHECK  FOR  CLOSURE  OF  NEW  ITERATED 
6  PARAMETER  AND  PREVIOUS 
3  PARAMETER  VALUE 


NO  CLOSURE 


CLOSURE 


STOP 


Figure  22.   Flow  Chart  of  First  Approach  to  Synthetic  Self-Calibration 
of  a  Gravity  Model:   Friction  Factor  B  Parameter  as  a 
Function  of  Areawide  Average  Travel  Time 
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gravity  models  are  shown  in  the  Tables  and  are  compared  to  those 
values  obtained  from  the  traditional  gravity  model  calibration  pro- 
cedure.  As  can  be  seen  in  Tables  25  and  26,  the  differences  in  B 
parameters  range  from  a  low  of  -0.0017  for  Anderson's  HBW  trips  to  a 
high  of  -0.1516  for  Anderson's  NHB  trips.   Lafayette's  largest  3 
parameter  difference  was  -0.0405  for  the  HBW  trips. 

The  total  internal  trip  tables  of  Lafayette  and  Anderson  as 
produced  by  the  test  gravity  model  were  then  assigned  to  their 
respective  networks  on  an  all-or-nothing  basis.   Volume  group  link 
comparisons  were  then  performed  with  their  corresponding  origin- 
destination  survey  all-or-nothing  traffic  assignments.   Tables  27  and 
28  present  the  results  of  these  volume  group  link  comparisons  for 
Lafayette  and  Anderson  respectively.   Total  link  trip  percent  root- 
mean-square-error  (PRMSE)  for  Lafayette  was  21.09  percent  as  compared 
to  the  12.67  percent  obtained  in  the  comparison  of  the  origin- 
destination  survey  assignment  with  the  traditionally  calibrated 
gravity  model  procedure.   The  corresponding  total  link  trip  PRMSE 
value  for  Anderson  was  28.15  percent  for  the  test  gravity  model  as 
compared  to  13.28  percent  for  the  traditionally  calibrated  model. 
Figures  23  and  24  graphically  portray  these  comparisons,  respectively 
for  Lafayette  and  Anderson. 

As  a  final  check  of  the  accuracy  of  the  test  gravity  model,  a 
link  by  link  analysis  of  the  model/0-D  link  load  ratios  was  performed 
for  both  Lafayette  and  Anderson.   For  Lafayette,  no  link  was  found  to 
have  an  unacceptable  link  load  ratio  while  Anderson  had  29  links,  or 
1.9  percent,  with  unacceptable  ratios. 
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Table  27.   Volume  Group  Link  Comparison  for  Lafayette  All-or-Nothing 
Traffic  Assignment:   Areawide  Self-Calibrating  Gravity 
Model  Versus  Origin-Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRl^lSE 

0-   500 

216 

12.20 

31.63 

47.61 

500-  1000 

116 

57.09 

20.25 

20.71 

1000-  2000 

197 

118.64 

14.43 

16.70 

2000-  3000 

140 

241.14 

11.25 

15.42 

3000-  5000 

172 

419.17 

8.91 

16.33 

5000-10000 

163 

719.53 

6.54 

15.33 

10000-15000 

31 

1409.39 

5.13 

14.23 

15000-20000 

1 

978.00 

4.35 

5.63 

Total 

1036 

290.00 

— — 

21.09 

Excludes  centroid  links 


Table  28.   Volume  Group  Link  Comparison  for  Anderson  All-or-Nothing 
Traffic  Assignment:   Areawide  Self-Calibrating  Gravity 
Model  Versus  Origin-Destination  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

604 

23.07 

34.63 

47.07 

500-  1000 

211 

116.38 

20.25 

25.15 

1000-  2000 

228 

268.23 

14.43 

23.60 

2000-  3000 

169 

345.47 

11.25 

18.39 

3000-  5000 

185 

571.16 

8.94 

18.28 

5000-10000 

108 

998.68 

6.54 

18.26 

10000-15000 

2 

1028.50 

5.12 

11.68 

Total 

1507 

247.34 

— 

28.15 

Excludes  centroid  links 
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Based  on  these  two  test  cases,  it  may  be  concluded  that  the 
self-calibrating  B  parameter  procedure  using  Equations  15,  16  and  17, 
is   not  as  reliable  a  procedure  as  the  traditionally  calibrated 
gravity  model  technique.   As  evidenced  by  the  Anderson  results,  the 
error  produced  by  this  approach  to  synthetic  self-calibration  can  be 
relatively  large  when  compared  to  the  traditional  calibration  results. 

Some  of  the  possible  reasons  for  the  poor  performance  of  the 
procedure  discussed  above  for  areawide  self-calibration  are  listed 
below. 

1.  Spatial  distribution  of  the  land  use  activities  and  the 
characteristics  of  the  highway  network  can  differ  significantly 
among  the  small  urban  areas.   Therefore,  the  factors  and 
decision  variables  considered  by  the  trip  maker  in  basing  his 
choice  of  destination  can  also  differ  with  respect  to  location. 

2.  Average  travel  time  for  an  entire  urban  area  is  a  highly 
aggregated  statistic  and  may  not  reflect  adequately  the 
transportation  system  options  encountered  by  the  individual 
trip  makers  at  various  locations  within  the  urban  area. 

3.  Terminal  times,  as  shown  in  Chapter  III,  can  vary  greatly 
among  the  study  areas.   This  variability  in  terminal  times 
is  directly  reflected  in  the  average  travel  times  for  the 
urban  areas.   This  complicates  the  regression  procedure  from 
which  the  calibration  equations  were  derived  by  adding 
additional  subjective  variance  to  the  statistical  curve 
fitting  process. 

The  above  three  items  were  considered  in  the  development  of  the  second 
approach  to  synthetic  self-calibration  of  a  gravity  model. 
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Orlsln  Zone  Specific  Approach  to  Synthetic  Self-Calibration 

of  a  Gravity  Model 
Table  29  compares  the  areawide  travel  time  statistics  for 
Lafayette  and  Anderson.   The  Table  shows  that  the  average  travel  time 
values  by  trip  purpose  are  approximately  equal  but  that  the  travel 
time  variances  are  completely  different  for  the  two  urban  areas.   The 
Anderson  travel  time  variances  are  approximately  one  half  of  those  of 
Lafayette.   However,  the  regression  analysis  used  to  develop  the 
areawide  B  parameter  calibration  equations  in  the  areawide  approach  did 
not  show  travel  time  variance  to  be  a  significant  factor  in  explaining 
any  of  the  variance  between  the  B  parameters  computed  on  the  basis  of 
tlie  data  from  eighteen  urban  areas.   The  analysis  conducted  for  the  de- 
velopment of  the  areawide  approach  and  the  possible  reasons  given  for 
the  poor  performance  of  the  derived  calibration  equations   suggest 
the  possibility  that  the  variance  of  data  within  an  urban  area  may  be 
more  significant  to  gravity  model  self-calibration  than  the  variance 
between  urban  areas. 


Table  29.   Lafayette  and  Anderson  Areawide  Travel  Time  Statistics 
Comparison 


Trip  Purpose 

HBW 
HBO 
NHB 


La 

tfayi 

ette 

Ande 

rson 

Mean 

Variance 

Mean 

Va 

riance 

10.60 

23.88 

11.44 

9.61 

8.57 

19.90 

10.04 

9.53 

9.17 

15.07 

9.50 

7.81 
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In  order  to  reveal  the  extent  to  which  the  travel  time  statistics 
may  vary  within  an  urban  area.  Table  30  was  constructed  using  the 
data  sets  of  Lafayette  and  Anderson.   Origin  zone  specific  average 
travel  times  and  travel  time  variances  were  computed  from  gravity 
models  with  constant  friction  factors,  that  is  the  opportunity  trip 
distribution  as  defined  by  Voorhees  (19).   As  presented  in  the  Table, 
origin  zone  average  travel  times  and  variances  possess  great  variability 
within  an  urban  area.   In  order  to  provide  a  better  pictorial 
representation  of  the  relationship  between  origin  zone  opportunity 
average  travel  time  and  zone  location,  Figures  25,  26  and  27  were  pre- 
pared for  the  HBW,  HBO,  and  NHB  trips  of  Lafayette,  Indiana.   These 
Figures  show  that  destination  opportunities  vary  with  both  location 
and  trip  purpose,  as  represented  by  the  origin  zone  opportunity  average 
travel  times,  and  therefore  it  is  reasonable  to  assume  that  the  desti- 
nation choice  process  itself  is  a  function  of  origin  zone  location. 


Table  30.   Origin  Zone  Opportunity  Travel  Time  Statistics  for 
Lafayette  and  Anderson 


Trip  Purpose 


Lafayette 


Anderson 


HBW 


Mean 


Min. 


Max. 


Average 

12.806 

8.759 

22.024 

13.075 

9.079 

22.016 

Variance 

14.497 

7.388 

22.755 

5.445 

2.475 

9.903 

HBO 

Mean 

Min. 

Max. 

Mean 

Min. 

Max. 

Average 

12.482 

8.520 

22.014 

13.160 

10.068 

21.198 

Variance 

15.137 

5.207 

22.693 

8.426 

4.603 

13.198 

NHB 

Mean 

Min. 

Max. 

Mean 

Min. 

Max. 

Average 

12.442 

8.358 

21.744 

12.930 

10.060 

20.815 

Variance 

13.910 

5.358 

21.257 

8.246 

4.453 

25.549 

Mean 


Min. 


Max. 
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Figure  25.   Origin  Zone  Opportunity  Average  Travel  Time  Values  for 
Lafayette  Home  Based  Work  Trips 
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Figure  26.   Origin  Zone  Opportunity  Average  Travel  Time  Values  for 
Lafayette  Home  Based  Other  Trips 
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Figure  27.   Origin  Zone  Opportunity  Average  Travel  Time  Values  for 
Lafayette  Non-Home  Based  Trips 
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What  is  being  proposed  is  a  stratified  approach  to  trip  distribu- 
tion.  Trip  generation,  the  modelling  process  performed  before  trip 
distribution,  has  evolved  from  zonal  regression  equations  to 
the  disaggregated  approach  of  cross  classification  of  residential  and 
non-residential  trip  rates.  Not  only  does  this  technique  simplify  data 
collection  and  improve  prediction  results,  but  it  has  also  been  found 
that  these  trip  rates  can  be  transferred  more  easily  from  one  study 
area  to  another  (10).  The  same  approach  may  be  followed  for  friction 
factors,  that  is,  friction  factors  or  friction  factor  equation 
parameters  may  be  considered  on  a  disaggregated  basis  by  computing 
these  values  by  zone  rather  than  using  an  aggregated  value  for  the 
entire  study  area. 

In  order  to  develop  a  disaggregated  synthetic  gravity  model,  it 
is  first  necessary  to  generate  the  data  on  which  this  development  is 
to  be  based.   The  study  areas  of  Lafayette,  Anderson,  and  Muncie, 
Indiana  were  used  for  this  purpose  because  of  data  availability  and 
the  wide  sampling  space  provided  by  the  differing  characteristics  of 
these  areas.   To  generate  the  data  base,  the  gravity  models  of  these 
three  urban  areas  were  calibrated  at  the  zonal  level.   This  was 
accomplished  by  using  the  origin-destination  survey  to  compute  origin 
zone  average  travel  times  and  then  using  these  values  to  repeatedly 
adjust  the  origin  zone  specific  parameters  of  the  negative  exponential 
friction  factor  function  until  convergence  of  the  gravity  model 
origin  zone  mean  travel  times  occurred.   The  data  generated  by  this 
procedure  consisted  of  the  following  for  each  origin  zone  and  for  each 
trip  purpose  -  HBW,  HBO,  and  NHB  trips. 
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3.    =  negative  exponential  friction  factor  parameter  for  zone  i 

of  the  origin  zone  calibrated  gravity  model, 
0.    =  trips  originating  from  zone  i, 
AVEO .  =  opportunity  average  travel  time  for  origin  zone  i  when 

3.  equals  0.0,  that  is,  friction  factor  of  constant 

value  1.0, 
VARO.  =  opportunity  travel  time  variance  for  origin  zone  i  when 

B-  equals  0.0, 
AVE.   =  origin  specific  average  travel  time  for  zone  i  of 

calibrated  gravity  model  (same  as  origin-destination 

survey  value) , 
VAR.   =  origin  specific  travel  time  variance  for  zone  i  of 

calibrated  gravity  model,  and 
TERM.  =  origin  zone  terminal  time. 

The  above  variables  were  then  used  in  a  curve  fitting  regression 
procedure  so  as  to  find  a  relationship  of  the  following  form: 


3.  =  f(AVEO.,  AVE.,  VARO.,  VAR.,  TER>I.) 
1         1      1       11       1 

Regression  Analysis  for  Development  of  an  Origin  Zone 
Specific  3  Parameter  Equation 
Having  generated  the  data  base  for  the  equation  development,  a 
weighting  scheme  was  employed  to  the  data  before  the  curve  fitting 
regression  analysis  was  begun.   Weighting  of  the  data  base  was  per- 
formed at  two  different  levels.   First,  the  data  was  weighted  by  the 
variable  0.,  that  is,  the  origin  zone  productions  (or  origins). 
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Neglecting  to  apply  this  weighting  factor  would  have  resulted  in 
zones  \;ith  only  a  few  trip  productions  (origins)  to  have  the  same 
effect  in  the  curve  fitting  regression  analysis  as  a  zone  of  a 
thousand  trips  or  more.   Secondly,  because  the  total  number  of  trips 
by  purpose  are  not  the  same  for  each  urban  area,  a  weighting  scheme 
was  applied  to  balance  this  inequality  and  thus  bias  toward  any  one 
particular  urban  area  was  minimized. 

As  was  noted  in  Chapter  III,  terminal  times  can  vary  greatly 
between  urban  areas.   Because  this  value  is  assigned  primarily  by 
subjective  judgment  and  because  origin  zone  terminal  times  must  be 
experienced  by  a  trip  maker   irrespective  of  the  location  of  the 
destination  choice,  it  was  decided  to  remove  origin  zone 
terminal  time   from  the  origin  zone  average  travel  time.   Origin  zone 
terminal  time,   however,  was   allowed  to  be  a  separate  variable  in 
the  regression  procedure. 

Finally,  dummy  variables  representing  the  three  urban  areas  were 
added.   This  allows  any  unique  and  significant  combinatorial  effects 
of  unspecified  variables  of  an  urban  area  to  be  represented  in  the 
analysis.   Thus  dummy  variables  make  major  differences  between  urban 
areas  to  be  apparent  in  the  analysis. 

The  initial  regression  analysis  resulted  in  an  equation  form 
possessing  very  acceptable  statistics  for  the  curve  fitting  of  the 
origin  zone  specific  6  parameters  to  their  respective  zonal  travel 
time  statistics  of  mean  and  variance.   However,  inspection  of  the 
residual  plot  revealed  a  greater  than  expected  number  of  outliers. 
Careful  review  of  the  zonal  data  and  the  previously  presented  sensitivity 


89 


analysis,  resulted  in  the  removal  of  the  variable  VAR^  from  the 
analysis,  that  is,  the  origin  zone  travel  time  variance.   A  close 
examination  of  the  data  revealed  that  zonal  travel  time  variance  does 
not  always  change  monotonically  with  a  monotonically  changing  friction 
factor  3  parameter.   For  some  zones,  travel  time  variance  first  in- 
creased and  then  decreased  for  an  increasing  B  parameter,  that  is, 
from  6  equal  to  zero  (opportunity  distribution)  and  upwards.   This 
phenomenon  is  probably  caused  by  the  gravity  model  attraction  constraints 
and  its  iteration  process. 

Having  removed  the  variable  VAR.,  further  regression  analysis 
resulted  in  an  equation  form  with  acceptable  statistics  and  residual 
plots.   However,  the  incorporation  of  the  3  parameter  equation  into  a 
gravity  model  did  not  provide  acceptable  results.   The  6  parameters 
were  found  to  be  converging  to  very  erroneous  values  with  many   not 
converging  at  all.   A  review  of  the  data  revealed  that  the  regression 
equations  developed  were  not  equations  of  how  the  calibrated  origin 
zone  B  parameters  related  to  various  variables,  but  represented 
actually  a  generalized  response  surface  of  any  B  parameter  with 
respect  to  the  independent  variables.   For  example,  Hyman  (9)  showed 
that  average  travel  time  varies  linearly  with  the  B  parameter,  and 
therefore  the  calibrated  origin  zone  6  parameter  is  just  one  point  on 
the  line.   To  eliminate  the  difficulty  with  a  generalized  response 
surface,  an  additional  restriction  v^7as  imposed  eliminating  the  in- 
clusion of  both  AVE.  and  AVEO .  in  the  regression  equation  simultaneously. 

11 
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This  additional  restriction  greatly  simplified  the  number  of 
possible  variable  combinations  to  be  studied.   Given  these  restrictions, 
the  host  regression  equations  that  could  be  developed  had  R-squared 
values  of  0.14,  0.52,  and  0.62  for  the  HBW,  HBO,  and  NHB  trips 
respectively.   Some  '-significance  was  also  indicated  by  the  urban  area 
dummy  variables;  a  result  totally  unwanted  for  it  is  necessary  to  have  a 
technique  applicable  to  all  urban  areas.   Because  of  the  low  R-squared 
values,  especially  for  the  HBW  trip  purpose,  and  tlie  significance  of 
the  dummy  variables,  these  equations  were  judged  unacceptable. 

It  had  been  initially  thought  that  the  development  of  origin  zone 

specific  self-calibrating  friction  factor  6  parameter  equations  was 

possible.   However,  for  reasons  cited  above,  this  approach  to  zonal 

calibration  was  abandoned.   The  results  of  the  regression  analysis  did, 

however,  suggest  another  possible  method  for  zonal  self-calibration. 

The  correlation  matrix  generated  by  the  regression  procedure  showed 

that  a  liigh  correlation  existed  between  the  variables  AVEO .  and  AVE.. 
^  11 

It  is  this  high  correlation  on  which  a  revised  second  approach  to 
synthetic  self-calibration  of  a  gravity  model  was  based. 

Revised  Second  Approach  to  Synthetic  Self-Calibration 

of  a  Gravity  Model 

The  initial  zone  specific  approach  to  estimate  the  S  parameter  of 

tlie  negative  exponential  friction  factor  equation  using  a  direct 

iterative  technique  proved  unsuccessful.   However,  a  high  correlation 

observed  between  AVEO.  and  AVE.  indicated  that  an  indirect  technique 

1        1 

to  estimating  the  ;-  parameter  could  be  possible.   Consequently,  a  re- 
vised procedure  was  developed  with  AVEO.  as  the  key  variable. 
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As  previously  mentioned,  AVEO.  is  the  opportunity  average  travel 

time  for  an  origin  zone,  as  defined  by  Voorhees  (19)  ,  minus  terminal 

time  for  the  origin  zone.   Origin  zone  terminal  time  was  subtracted  so 

as  to  separate  it  from  the  link  travel  time  and  destination  zone 

terminal  time  which  are  themselves  a  function  of  destination  zone 

location.   The  high  correlation  of  AVE.  with  AVEO.  suggests  that  the 

destination  choice  process  is  closely  related  to  the  geographic 

location  of  the  urban  area  activity  system  with  respect  to  the  origin 

zone.   It  was  felt  that  if  an  estimating  equation  for  AVE.  could  be 

developed  as  a  function  of  AVEO.  and  TEPJI.,  then  these  equations  could 

11 

be  incorporated  in  a  gravity  model  for  synthetic  self-calibration. 

Regression  analysis  of  AVE.  on  AVEO.  and  TERM,  was  accomplished 
by  weighting  the  zonal  data,  as  previously  described,  in  order  to 
properly  balance  the  data  by  trip  productions  (origins)  and  by  trip 
purpose.   Equations  18,  19,  and  20  present  the  regression  analysis 
results  for  the  HBW,  HBO,  and  MHB  trips  respectively.   Total  R-squared 
and  the  change  in  R-squared  for  the  independent  variables  are  also 
given. 


HBW:   AVE.  =  -0.4328+1.2761  TERI-1.+0. 8921  (AVEO. -TERM. ) 
1  111 

R^  =  .865        AR*^  =  .114       AR^  =  0.751      (18) 


HBO:   AVE.  =  -0.7499+1.3548  TERl-l.+O.  7558  (AVE0.-TER.M. ) 
1  111 

R^  =  .848        AR^  =  .187       AR^  =  .661       (19) 


NHB:   AVE.  =  0.5245+1.3416  TERM. +0.6073  (AVEO .-TERM. ) 
1  111 

R^  =  .719        AR^  =  .423       AR^  =  .296       (20) 
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Thf  L'stiniating  equations  possess  R-squared  values  higher  than 
those  achieved  in  earlier  attempts.   The  strong  statistical  fit  of 
the  equations  clearly  shows  the  importance  of  AVEO.,  the  opportunity 
average  travel  time  on  the  destinati(5n  choice  process.   This  in  turn, 
demonstrates  the  large  effect  that  origin  zone  location  and  the 
distribution  of  land  use  activities  about  the  origin  zone  have  on  the 
destination  choice  process. 

Relative  importance  of  travel  time  by  trip  purpose  in  the 
destination  choice  process  is  indicated  by  the  coefficients  of  the 
AVEO.  variable.   In  general,  the  home  based  work  (HBW)  trips  do  not 
have  much  choice  in  terms  of  destination,  the  home  based  other  (HBO) 
trips  have  some  or  many  choices  of  destination,  while  the  non-home 
based  (NHB)  trips  have  many  choices  of  destination.   It  can  be  seen 
from  the  Equations  that  the  importance  of  travel  time  is  the  lowest 
for  HBW,  highest  for  NHB,  and  of  some  intermediate  magnitude  for  HBO 
trips.   This  relative  importance  is  also  reflected  in  the  corresponding 
change  of  the  R-squared  value  for  the  AVEO.  variable. 

In  addition,  variable  TERM.,  the  origin  zone  terminal  time,  be- 
comes a  greater  proportion  of  total  travel  time  as  trip  purpose 
changes  from  HBW,  to  HBO,  and  to  NHB,  as  is  evidenced  by  the  change  in 
R-squared.   Tliis  is  because  terminal  times  are  zone  specific  constants 
and  cannot  be  reduced  by  destination  choice  unlike  the  link  travel 
t  imes . 

Also  of  interest  are  the  coefficients  of  the  origin  zone  terminal 
times.   It  was  first  expected  that  the  coefficients  would  have  a 
magnitude  very  close  to  one.   The  fact  that  the  coefficients  are 
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greater  than  i.Omay  possibly  be  attributable  to  that  constant  part  of 
the  destination  zone  terminal  time  common  to  all  zones,  that  is,  the 
minimum  terminal  time. 

Having  developed  estimating  equations  by  trip  purpose  for  the 
origin  zone  average  travel  time,  AVE.,  a  procedure  was  then  formulated 
to  incorporate  these  equations  into  a  gravity  model. 

Incorporation  of  the  Estimating  Equations  Into  a  Gravity  Model 
One  method  of  calibrating  the  gravity  model  would  be  to  first 
distribute  all  trips  with  a  constant  friction  factor,  thereby  allowing 
for  the  computation  of  AVEO.,  the  origin  zonal  specific  opportunity 
average  travel  time.   Origin  zone  average  travel  time,  AVE.,  could 
then  be  computed  using  the  regression  equations  discussed  in  the 
previous  section.   These  origin  zone  average  travel  times  then  could 
be  used  to  repeatedly  adjust  the  origin  zone  specific  6  parameters 
until  the  gravity  model  origin  zone  average  travel  times  match  the 
estimated  values.   However,  this  type  of  calibration  procedure  would 
require  many  iterations  before  convergence  for  all  the  zones  occur. 
The  many  iterations  required  for  convergence  would  consume  a  large 
amount  of  computer  time.   This  would  be  an  undesirable  characteristic 
for  a  synthetic  model  to  possess. 

Consequently,  another  method  was  developed  for  eventual  adoption 
for  synthetic  calibration  and  was  based  on  three  assumptions  regarding 
the  interaction  of  origin  zone  average  travel  time  and  the  zonal  6 
parameter.   The  first  assumption  was  that  the  origin  zone  average 
travel  time  varies  linearly  with  the  origin  zone  p.  parameter  (9). 
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Second,  zonal  average  travel  time  is  independent  of  the  p  parameters 
of  otiier  zones.   Third,  zonal  average  travel  time  is  independent  of 
the  attraction  constraint  iteration  procedure.   The  key  assumption  of 
linearity  does  approximate  reality  as  can  be  evidenced  in  Figures  28 
and  29  by  the  plots  of  origin  zone  average  travel  time  versus  areawide 
B  parameter,  that  is  all  zonal  3  parameters  equal.   Figure  28  contains 
the  plots  of  zones  41  through  49  for  the  HBVJ  trips  of  Lafayette.   The 
Figure  shows  typical  plots  in  which  the  zonal  average  travel  time 
approximately  varies  linearly  with  the  3  parameter.   Figure  29  con- 
tains plots  of  zones  91  through  99,  also  for  Lafayette  HBW  trips. 
This  Figure  shows  greater  nonlinearity  for  the  Lafayette  HBW  trips 
but  even  these  are  not  great  enough  to  reject  the  general  assumption 
of  linearity,  especially  when  one  considers  that  most  ,'  parameters  are 
less  than  0.40  and  that  zonal  results  are  aggregated  up  to  the  study 
area  level.   These  assumptions,  however,  do  not  always  hold  true  for 
some  of  the  individual  zones  when  the  origin  zone  >  parameters  are 
allowed  to  vary  independently  as  in  the  proposed  model.   This  aspect 
could  not  be  detected  until  after  the  model  was  developed  and  tested, 
but  aggregation  of  the  results  up  to  the  entire  study  area  nullified 
these  minor  errors. 

The  estimating  equations  of  AVE.,  zonal  average  travel  time,  and 
ttie  three  assumptions  listed  above  were  incorporated  into  a  gravity 
model  procedure  that  can  best  be  described  in  conjunction  with  Figure 
30.   First,  trips  are  distributed  with  all  zonal  ?    parameters  equal  to 
0.0  (constant  friction  factor)  and  then  the  resulting  zonal  opportunity 
average  travel  times,  AVEO.  are  computed.   Second,  trips  are  distributed 
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STEP  I 


0.0 


^,-      SET  VALUE 
Origin  Zone  3  Parameter 


Step  1  -  Average  at  6  equal  zero 
(0.0)  -  AVEO. 


Average  at  6  equal  to  a 
set  value. 


Step  2 

Step  3  -  Average  from  regression 
equation,  AVE. 

Step  4  -  Interpolate  for  new  3 


Figure  30.   Generalized  Flow  Chart  of  Synthetically  Self-Calibrated, 

Origin  Zone  Specific  Gravity  Model  as  Applied  to  a  Typical 
Zone 


with  all  origin  zone  3  parameters  equal  to  a  preselected  value  and  then 
the  origin  zone  average  travel  times  are  computed.   Third,  using  the 
estimating  equation  for  the  appropriate  trip  purpose  and  the  AVEO. 
values  computed  in  step  one,  the  origin  zone  average  travel  times, 
AVE.,  are  determined.   Fourth,  using  the  data  points  from  steps  one 
and  three  and  the  estimated  origin  zone  average  travel  times,  AVE., 
the  zonal  3  parameters  are  obtained  through  interpolation.   All  trips 
are  then  distributed  using  the  interpolated  origin  zone  specific  6 
parameters. 
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Evaluation  of  the  Proposed  Model 
The  gravity  model  procedure,  as  just  described,  was  used  to 
distribute  the  trips  for  botli  Lafayette  and  Anderson.   The  resulting 
average  travel  time  statistics  by  trip  purpose  and  for  all  internal 
trips  are  given  in  Tables  31  and  32  for  Lafayette  and  Anderson 
respectively.   Comparison  of  these  statistics  to  the  corresponding 
origin- destination  survey  values  is  also  shown.   Lafayette's  total 
internal  trip  average  travel  time  was  only  0.71  percent  greater  than 
tlie  0-D  survey  value,  while  Anderson's  total  internal  average  travel 
time  was  identical  to  the  0-D  survey  value. 

The  average  travel  times  by  trip  purpose  also  indicated  very 
acceptable  results;  the  largest  percentage  error  for  either  urban 
area's  individual  trip  purposes  was  3.78  percent,  this  being  the  only 
one  outside  of  FHWA's  suggested  +3  percent  range. 

Volume  group  link  comparisons  of  the  two  synthetic  gravity  models' 
all-or-nothing  traffic  assignments  were  performed  using  their 
corresponding  0-D  survey  traffic  assignments  as  the  base  for  comparison. 
The  PRMSE  values  by  volume-group  are  presented  in  Tables  33  and  34  for 
Lafayette  and  Anderson  respectively.   Total  link  trip  PRMSE  for 
Lafayette  was  computed  to  be  12.67  percent.   This  value  is  practically 
identical  to  the  corresponding  value  (12.48)  obtained  from  the  gravity 
model  calibrated  by  the  traditional  procedure.   Total  link  trip  PRIISE 
was  calculated  to  be  14.55  percent  as  compared  to  the  traditionally 
calibrated  gravity  model  value  of  13.28  percent  for  Anderson.   Volume 
group  link  comparison  statistics  of  Tables  33  and  34  are  also  shown 
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Table  31.   Origin  Zone  Specific,  Synthetically  Self-Calibrating  Gravity 
Model  Travel  Time  Statistics  for  Lafayette 


Trip  Purpose 


HBW 
HBO 
NHB 
Total 


Average  Travel  Time 
Survey        Model 


0.57 

10.97 

8.56 

8.65 

9.12 

8.93 

9.05 

9.09 

Percent  Difference 


+  3.78 
+1.05 
-2.08 
+0.71 


Table  32.   Origin  Zone  Specific  Synthetically  Self-Calibrating  Gravity 
Model  Travel  Time  Statistics  for  Anderson 


Trip  Purpose 


Average  Travel  Time 
Survey        Model 


Percent  Difference 


HBW 
HBO 
NHB 
Total 


11.44 

11.27 

10.04 

10.02 

9.50 

9.64 

10.14 

10.14 

-1.49 

-0.20 

+1.47 

0.00 
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Table  33.   Lafayette  All-or-Nothing  Traffic  Assigninent  Comparison: 

Origin  Zone  Specific  Synthetically  Self-Calibrating  Gravity 
Model  Versus  0-D  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PR14SE 

PRMSE 

0-   500 

216 

11.59 

34.63 

50.15 

500-  1000 

116 

44.09 

20.25 

18.58 

1000-  2000 

197 

38.88 

14.43 

14.56 

2000-  3000 

140 

81.90 

11.25 

11.20 

3000-  5000 

172 

177.01 

8.94 

11.46 

5000-10000 

163 

209.29 

6.54 

8.78 

10000-15000 

31 

444.58 

5.12 

6.29 

15000-20000 

1 

-9.00 

4.35 

0.05 

Total 

1036 

101.43 

— 

12.67 

Table  34.   Anderson  All-or-Nothing  Traffic  Assignment  Comparison: 

Origin  Zone  Specific  Synthetically  Self-Calibrating  Gravity 
Model  Versus  0-D  Survey 


No.  of 

Mean 

Expected 

Observed 

Volume  Group 

Links 

Difference 

PRMSE 

PRMSE 

0-   500 

604 

16.15 

34.63 

42.27 

500-  1000 

211 

49.48 

20.25 

17.94 

1000-  2000 

228 

47.75 

14.43 

12.72 

2000-  3000 

169 

-43.60 

11.25 

9.82 

3000-  5000 

185 

-115.35 

8.94 

9.13 

5000-10000 

108 

-274.01 

6.54 

8.16 

10000-15000 

2 

-1201.50 

5.12 

13.71 

Total 

1507 

-17.68 

_  — 

14.55 
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graphically  in  Figures  31  and  32.   Also  shown  on  tlie  Figures  are  tlie 
expected  PRMSE  values  associated  with  a  ten  percent  dwelling  unit 
sampling  rate  for  the  0-D  home  interview  survey.   The  PRMSE  values 
for  both  Lafayette  and  Anderson  were  observed  to  be  very  acceptable. 

As  a  final  test  of  the  developed  procedure  for  synthetic  self- 
calibration  of  a  gravity  model,  a  model/0-D  link  load  ratio  analysis 
was  performed  so  as  to  determine  the  number  of  links  outside  of  the 
suggested  ratio  ranges.   For  the  Lafayette  urban  area,  no  model/0-D 
link  load  ratio  was  observed  to  fall  outside  of  the  suggested  ranges. 
For  the  Anderson  study  area  only  two  links  indicated  unacceptable 
model/0-D  link  load  ratios.  However,  these  two  links  are  the  same 
links  as  those  found  unacceptable  in  the  traditionally  calibrated 
gravity  model  traffic  assignment  and  are  of  relatively  low  volumes. 

Because  of  the  excellent  results  obtained  from  this  approach  for 
gravity  model  calibration  on  the  Lafayette  and  Anderson  study  areas, 
the  program  logic  was  rewritten  so  as  to  incorporate  it  into  program 
UMODEL  of  the  UTPS  package  for  additional  testing  on  the  Indiana  State 
Highway  Commission's  IBM  370  computer.   Appendix  C  contains  the  listing 
of  this  program  code.   This  allowed  testing  of  the  program's  modelling 
capability  on  the  Muncie,  Indiana  data  set  and  provided  for  the 
possible  quick  implementation  of  the  proposed  gravity  model  into  the 
normal  planning  process  if  these  tests  proved  successful. 

Average  travel  time  statistics  by  trip  purpose  as  computed  from 
the  running  of  the  test  gravity  model  are  given  in  Tabic  35  for  the 
Muncie  urban  area.   Also  shown  in  the  Table  are  the  0-D  survey  average 
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travel  times  and  the  percent  differences.   Total  trip  average  travel 
time  as  computed  by  the  test  gravity  model  underestimated  the  0-D  survey 
value  by  only  a  -0.80  percent,  with  the  largest  trip  purpose  error 
being  a  -2.79  percent  for  HBW  trips.   The  errors  are  far  within 
acceptable  limits. 


Table  35.   Travel  Time  Statistics,  Generated  by  Origin  Zone  Specific 
Synthetic  Gravity  Model  for  the  Muncie  Urban  Area 


Areawide  Average  Travel  Time  „ 
° ■ Percent 

Trip  Purpose        0-D  Survey  Model  Difference 


HBW  9.07  8.82  -2.79 

HBO  7.48  7.37  -1.46 

NHB  6.93  7.04  +1.63 

Total  7.52  7.46  -0.80 

The  internal  trip  tables  as  generated  by  the  test  gravity  model 
were  next  assigned  to  Muncie's  network  on  an  all-or-nothing  basis. 
Internal  trip  tables  from  the  0-D  survey  were  also  assigned  to  the 
Muncie  network  on  an  all-or-nothing  basis  so  that  a  comparison  between 
the  model  and  0-D  survey  assignments  could  be  performed. 

The  resulting  PRMSE  values  from  a  volume  group  link  comparison  of 
these  tvs^o  assignments  were  found  to  be  much  higher  than  expected.   It 
was  therefore  decided  to  further  evaluate  the  performance  of  the  model. 
For  this  purpose  a  gravity  model  calibration  for  the  HBW,  HBO,  and  NHB 
trips  was  undertaken  using  the  traditional  calibration  procedure,  with 
the  resulting  trip  tables  assigned  to  the  Muncie  network.   A  volume-group 
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link  comparison  of  Lhis  assignment  with  the  0-D  survey  assignment  was 
subsequently  performed.   In  this  manner,  the  modelling  capabilities 
of  the  origin  zone  specific  calibration  procedure  can  be  compared  with 
the  traditionally  calibrated  gravity  model  as  the  "standard"  model. 
If  the  proposed  procedure  would  model  trip  distribution  as  accurately 
as  the  traditional  model,  then  it  could  be  concluded  that  the  proposed 
procedure  is  probably  a  valid  synthetic  trip  distribution  model  for 
Muncie. 

The  volume-group  link  comparison  for  the  traditionally  calibrated 
gravity  model  resulted  in  a  total  link  trip  PRMSE  value  much  larger 
than  the  value  obtained  from  the  proposed  calibration  procedure.   A 
second  check  was  also  made  by  performing  the  same  analysis  on  the  trip 
tables  generated  by  the  Muncie  transportation  study,  which  also  in- 
cluded socio-economic  "K"  factors.   Total  link  trip  PRMSE  for  the  "K" 
factor  gravity  model  was  observed  to  be  slightly  higher  than  the 
traditionally  calibrated  gravity  model.   Table  36  contains  the  volume- 
group  PRMSE  values  for  all  three  gravity  models  for  Muncie  as  described 
above.   As  is  evident  from  the  Table,  the  origin  zone  specific  synthetic 
calibration  procedure  performed  much  better  than  the  traditionally 
calibrated  gravity  model  with  or  without  "K"  factors.   A  plot  of  the 
volume  group  link  comparison  PRMSE  values  resulting  from  the  origin 
zone  specific  synthetic  self-calibration  gravity  model  and  those  from 
the  traditionally  calibrated  model  is  given  in  Figure  33.   This  Figure 
illustrates  that  the  proposed  procedure  provides  more  accurate  results 
for  the  higher  volume  links  of  Muncie  as  compared  to  the  lower  volume 
links.   As  the  higher  volume  links  are  critical  in  the  process  of 
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transportation  facility  design,  the  proposed  form  of  gravity  model 
appears  to  be  an  effective  procedure. 

A  possible  explanation  for  the  high  PRMSE  values  for  all  three 
Muncie  models  may  be  the  fact  that  the  Muncie  study  area  included  a 
relatively  high  number  of  zones  and  links  for  an  urban  area  of  its 
population  size.   The  zones  and  links  are  proportionally  much  greater 
than  either  Lafayette  or  Anderson,  suggesting  that  the  Muncie  zone  and 
network  system  were  probably  too  detailed  for  the  number  of  trips 
produced.   For  example,  the  Muncie  data  included  1249  out  of  2698  links 
with  an  average  daily  traffic  of  less  than  only  500  vehicles  per  day; 
peak  hour  traffic  on  these  links  would  probably  not  exceed  50-60 
vehicles  per  hour.   Muncie  had  271  internal  zones  as  compared  to 
Anderson's  169  internal  zones  and  Lafayette's  111  internal  zones.   A 
large  number  of  zones  in  comparison  to  the  number  of  trips  can  also 
cause  numerical  problems  when  rounding  techniques  are  used  to  integerize 
the  zonal  interchanges.   Because  of  production  and  attraction  con- 
straints, rounding  techniques  may  accumulate  trips  from  some  zones  until 
a  whole  integer  trip  is  obtained.   Muncie's  288613  internal  vehicle 
trips  must  be  distributed  about  271  by  271  cells  of  the  trip  tables, 
or  an  average  of  3.9  trips  per  cell,  which  means  there  would  be  hundreds 
of  cells  with  less  than  one  trip  after  the  gravity  model  is  applied. 

Chapter  Conclusions 
In  this  Chapter,  the  development  of  an  improved  procedure  for 
synthetic  calibration  of  gravity  models  has  been  discussed.   The 
proposed  synthetic  procedure  Involves  an  origin  zone  specific 
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self-calibrating  gravity  model  In  whicli  the  only  input  data  required 
are  the  zonal  productions  and  attractions,  a  zone-to-zone  travel  time 
matrix  (skim  tree),  and  the  origin  zone  terminal  times.   A  travel  time 
distribution  is  no  longer  needed  for  calibration,  thereby  eliminating 
the  necessity  for  an  0-D  survey  for  internal-internal  trips. 

Tests  conducted  on  three  separate  study  areas  indicated  that  the 
proposed  model  is  able  to  reproduce  trip  patterns  as  accurately  as  the 
traditionally  calibrated  gravity  model  procedure  which  is  calibrated 
on  0-D  survey  data.   In  comparison  to  other  methods  of  synthetic  trip 
distribution,  the  proposed  gravity  model  provided  results  with  a  much 
greater  degree  of  reliability.   The  improved  reliability  was  achieved 
by  calibrating  the  model  at  the  zonal  level  rather  than  at  the 
aggregate  level  of  the  entire  study  area.   This  approach  minimizes 
overall  modelling  error  by  aggregating  the  individual  zonal  results 
up  to  the  urban  area. 

Development  of  the  proposed  procedure  was  also  based  on  the 
consideration  that  the  trip  distribution  is  critically  dependent  on 
spatial  distribution  of  land  use  activities  about  the  origin  zone. 
This  consideration  was  incorporated  in  the  proposed  procedure  through 
the  explicit  measurement  of  the  origin  zone  specific  opportunity 
travel  time  distribution.   The  opportunity  distribution  was  represented 
in  the  model  by  the  opportunity  average  travel  time.   It  is  this  key 
variable  on  which  the  model  was  developed  and  upon  which  the  very 
acceptable  model  results  are  due. 


110 


CHAPTER  V 
SUMMARY  AND  CONCLUSIONS 

The  present  study,  a  part  of  an  overall  research  project  involving 
synthetic  travel  demand  modelling,  included  only  the  trip  distribution 
phase  of  the  total  modelling  process.   Synthetic  trip  distribution  is 
a  method  by  which  the  travel  time  impedance  function  of  the  trip 
distribution  model  can  be  established  with  little  or  no  data  obtained 
directly  from  an  extensive  home- interview  travel  survey  of  the  area 
under  study.   This  Chapter  summarizes  the  results  and  conclusions  of 
this  synthetic  trip  distribution  study,  and  makes  recommendations  as 
to  what  technique  should  be  employed  in  synthetic  trip  distribution 
modelling  for  small  urban  areas  in  Indiana.   Recommendations  for 
further  research  are  also  presented  at  the  end  of  this  Chapter. 

The  great  importance  of  travel  time  impedance  factors  in  a  trip 
distribution  model  for  small  urban  areas  was  clearly  illustrated  by 
the  values  of  average  travel  time  and  the  traffic  assignment  statistics 
generated  from  trip  distribution  models  without  such  time  dependent 
factors.   Errors  in  the  resulting  travel  time  values  and  traffic 
assignments  were  of  a  magnitude  so  large  as  to  make  any  of  the 
generated  output  practically  useless.   The  analysis  using  gravity  models 
with  constant  friction  factors  for  the  entire  range  of  time  values 
clearly  establis'ned  the  need  and  importance  of  friction  factors. 
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It  is  important  for  systems  modellers  to  know  how  sensitive  a 
particular  model  is  with  respect  to  the  variable  parameters  of  the 
model.   Once  the  model  sensitivity  is  determined,  approximate  error 
limits  on  the  variable  parameters  can  be  estimated  or  the  modelling 
technique  can  be  altered  to  achieve  an  acceptable  level  of  accuracy. 
The  same  is  true  for  trip  distribution  modelling,  especially  when 
synthetic  metliods  are  to  be  used.   An  analysis  conducted  in  the 
present  study  indicated  that  the  sensitivity  of  the  gravity  model  to 
friction  factor  errors  may  be  highly  dependent  upon  the  particular 
urban  area  being  modelled.   That  is,  the  topographical  barriers,  land 
use  activity  distribution,  and  network  characteristics  interact 
uniquely  for  each  urban  area,  therefore  resulting  in  gravity  models 
with  varying  sensitivity.   As  was  demonstrated  in  the  present  study, 
the  gravity  model  for  the  Lafayette  study  area  was  found  to  be  very 
sensitive  to  deviations  of  the  calibration  parameters  in  the  travel 
time  impedance  function  from  the  actual  values  obtained  from  the 
origin-destination  survey.   Average  travel  times  by  trip  purpose  and 
traffic  assignment  link  volumes  changed  significantly  with  relatively 
small  deviations  of  the  calibration  parameter.   This  was  not  the  case, 
however,  for  Anderson  where  average  travel  times  and  link  volumes  did 
not  change  as  rapidly  as  they  did  for  Lafayette  when  the  calibration 
parameters  were  varied.   In  other  words,  the  calibration  parameter 
accuracy  required  for  the  generation  of  acceptable  model  results  is  a 
function  of  the  particular  urban  area  being  studied  and  therefore, 
cannot  be  generalized  for  any  urban  area. 
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The  borrowing  of  friction  factors  from  a  similar  urban  area  has 
long  been  suggested  as  a  technique  for  synthetic  gravity  models.   The 
most  often  cited  criterion  for  the  selection  of  a  similar  urban  area 
from  which  to  borrow  factors  is  population  size.   Population  size  is 
also  a  criterion  by  which  standardized  friction  factors  or  parameters 
are  categorized.   Both  of  these  techniques  therefore  assume  that  a 
correlation  exists  between  population  size  and  the  friction  factors 
themselves,  and  that  the  variability  with  respect  to  population  size 
is  acceptable.   As  it  was  shown  in  Chapter  III,  a  correlation  does 
exist,  but  a  large  variability  with  respect  to  population  was  observed. 
This  large  variability  therefore  weakens  the  confidence  in  the  results 
of  a  model  in  which  the  friction  factors  are  borrowed  or  standardized 
values  used. 

As  an  example  of  the  problems  that  may  occur  when  friction  factors 
are  borrowed  from  a  similar  urban  area,  the  calibrated  3  parameters 
for  Lafayette  and  Anderson,  Indiana  were  exchanged  between  the  two 
study  areas  of  roughly  the  same  population  size.   Average  travel  times 
by  trip  purpose  were  shown  to  possess  a  relatively  large  error, 
greater  than  the  FHWA's  suggested  +3.0  error  limit.   However,  off- 
setting errors  in  terms  of  trip  purpose  average  travel  times  resulted 
in  acceptable  traffic  assignments.   In  this  situation,  chance  and  not 
modelling  accuracy  resulted  in  acceptable  results. 

The  use  of  standardized  friction  factors  or  parameters  has  also 
been  suggested  in  the  literature  as  another  possible  synthetic  technique 
for  trip  distribution  modelling.   Standardized  a  parameters,  as  used 
in  the  inverse  travel  time  impedance  function,  were  developed  by  the 
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Comsis  Corporation  and  categorized  by  trip  purpose  and  population 
size  (3).   These  standardized  a   parameters  for  synthetic  trip  distri- 
bution modelling  were  tested  using  the  two  study  areas  of  Lafayette 
and  Anderson,  Indiana.   For  Lafayette,  the  average  travel  times  by 
trip  purpose  were  within  the  +3  percent  ranges  as  suggested  by  FHWA 
(13)  and  the  percent  root-mean-square-error  values  by  volume  group 
were  also  acceptable.   However,  12  of  the  1036  non-centroid  links  had 
unacceptable  model/0-D  survey  link  load  ratios.   Statistics  for 
Anderson  showed  that  the  non  home-based  trip  average  travel  time  was 
overestimated  by  5.74  percent.   However,  due  to  the  underestimate  of 
the  home  based  work  trip  average  travel  time  which  helped  to  offset 
the  non-home  based  travel  time  error  and  also  because  of  the  previously 
discusst^'d  insensitivity  of  the  Anderson  gravity  model  and  traffic 
assignment  to  friction  factor  errors,  an  acceptable  assignment  of  the 
traffic  to  the  network  was  observed.   Based  on  these  to  tost  cases, 
one  cannot  definitely  conclude  whether  or  not  the  synthetic  technique 
of  using  standardized  a  parameters  is  an  acceptable  procedure.   How- 
over,  because  of  the  12  links  of  the  Lafayette  network  with  unacceptable 
link  load  ratios,  the  high  overestimate  of  non-home  based  trip  average 
travel  time,  and  the  general  conclusion  that  urban  area  characteristics 
can  differ  greatly  suggests  the  possibility  that  standardized  a 
parameters  may  not  alv/ays  be  reliable. 

A  common  approach  to  represent  the  travel  pattern  of  an  urban 
area  is  to  use  the  travel  time  statistics  of  average  travel  time  and 
travel  time  variance  as  obtained  either  by  analyzing  the  0-D  survey 
data  conducted  in  the  study  area  or  as  generated  by  a  trip  distribution 
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model.   These  travel  time  statistics,  as  determined  from  a  past  0-D 
survey  of  one  urban  area,  are  sometimes  used  to  possibly  validate  the 
synthetic  trip  distribution  model  results  of  another  similar  urban 
area  being  studied.   Choosing  a  similar  urban  area  from  which  travel 
time  statistics  could  be  compared  has  mainly  been  based  on  the  criterion 
of  population  size.   This  criterion  is  based  on  the  assumption  that  a 
strong  correlation  exists  between  average  travel  time  and  population 
size.   It  has  been  observed  in  this  present  study  that,  although  there 
exists  a  correlation  among   the  average  travel  times  and 
population  size,  the  variability  of  average  travel  time  with  respect 
to  population  size  is  large.   Consequently,  the  validation  of  a 
syntlietic  trip  distribution  model  by  comparing  only  the  average  travel 
times  is  not  a  valid  approach.   Furthermore,  the  estimation  of  average 
travel  time  based  on  population  size  alone  would  be  o f  little  beneficial 
use  due  to  this  variability  in  the  data. 

The  same  problem,  as  discussed  above,  exists  for  the  travel  time 
frequency  distributions.   Even  though  the  general  shapes  of  the  fre- 
quency distributions  may  appear  quite  similar,  even  a  small  variation 
or  shift  in  a  distribution  can  significantly  alter  average  travel  time 
and  therefore  also  the  accuracy  of  the  model.   Again,  comparison  of 
travel  time  frequency  distributions  would  not  be  beneficial  in 
verifying  synthetic  trip  distribution  model  results. 

Because  of  the  inability  of  the  synthetic  techniques  of  borrowing 
friction  factors  or  using  standardized  friction  factor  values  to  con- 
sistently and  reliably  model  trip  distribution  travel  patterns,  a  large 
effort  was  therefore  undertaken  in  the  development  and  testing  of  an 
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acceptable  procedure  for  synthetic  self-calibration  of  gravity  models. 
It  was  first  found  that  a  gravity  model  with  an  areawide  friction 
factor  equation  parameter  as  a  function  of  areawide  travel  time 
statistics  could  not  adequately  model  travel  patterns.   The  next 
attempt  was  to  develop  an  origin  zone  specific  equation,  but  no 
satisfactory  regression  equation  could  be   developed  to  predict  the 
P  calibration  parameter  of  the  negative  exponential  function  on  the 
basis  of  origin  zone  travel  time  statistics.  The  final  approach  involved 
the  use  of  regression  analysis  of  origin  zone  mean  travel  time, 
terminal  time,  and  origin  zone  opportunity  mean  travel  time,  that  is 
the  average  travel  time  when  friction  factor  values  in  a  gravity  model 
are  held  constant.   The  statistical  result  was  satisfactory  and  the 
associated  procedure  developed  was  tested  on  the  study  areas  of 
Lafayette,  Anderson,  and  Muncie.   Comparisons  of  average  travel  times, 
percent  root-mean-square-error  values  by  volume  group,  and  model/O-D 
survey  link  load  ratios  indicated  very  acceptable  results  for  all  three 
study  areas.   The  accuracy  of  the  trip  distribution  results  were 
found   to  be  at  least  as  accurate  as  those  provided  by  the  traditionally 
calibrated  gravity  model. 

The  results  of  the  present  research  project  indicated  that  the 
most  reliable  synthetic  trip  distribution  modelling  technique  in- 
vestigated is  the  origin  zone  specific  synthetically  self-calibrating 
gravity  model.   This  teclinique  is  recommended  for  implementation 
because  of  its  ability  to  distribute  trips  at  an  accuracy  level  very 
close  or  better  than  the  traditional  technique  of  calibration  with 
extensive  origin-destination  data.   As  compared  to  other  synthetic 
techniques,  this  new  approach  has  several  advantages.   The  most 
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important  feature  of  the  proposed  procedure  is  its  more  direct 
linkage  to  the  spatial  distribution  of  land  use  activities  within  an 
urban  area  as  quantified  by  the  origin  zone  specific  opportunity 
average  travel  times.   This  inclusion  of  the  origin  zone  specific 
measure  of  spatial  land  use  activities  allows  the  gravity  model  to 
be  sensitive,  for  any  origin  zone,  to  the  surrounding  land  uses,  in 
addition  to  various  network  configurations.   This  procedure  is  there- 
fore probably  more  sensitive  than  any  areawide  model  to  changes  that 
may  occur  in  the  study  area,  such  as  a  new  shopping  center,  a  new 
manufacturing  plant  opening  or  a  major  network  alteration.   Another 
advantage  is  its  ease  of  application.   No  collection  of  origin- 
destination  survey  data  is  needed  nor  is  it  necessary  to  find  a 
similar  urban  area  from  which  to  borrow  information  or  to  which  model 
validation  comparisons  are  to  be  made.   As  this  present  study  showed, 
the  criterion  for  the  selection  of  a  similar  study  area  need  not 
guarantee  the  finding  of  such  a  similar  urban  area.   For  example,  the 
urban  areas  of  Lafayette,  Anderson,  and  Muncie  have  a  population  of 
approximately  the  same  size  but  were  found  to  have  many  differing 
characteristics.   The  proposed  model,  however,  successfully  simulated 
the  travel  patterns  of  these  three  urban  areas. 

Recommendations  for  Further  Research 
To  support  the  validity  and  applicability  of  the  recommended 
origin  zone  specific  approach  for  synthetic  self-calibration  of  a 
gravity  model,  further  test  applications  should  be  performed  on 
additional  transportation  study  areas.   This  would  require  the  applica- 
tion of  the  model  to  various  areas  from  different  regions  of  the  country. 
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In  order  to  simplify  such  application,  the  proposed  procedure  has 
been  coded  so  as  to  be  incorporated  into  LTMTA's  ITMODEL  program  and 
also  to  be  compatible  with  FHWA's  PLANPAC  package. 

Because  of  the  simplified  technique  used  in  the  gravity  model, 
it  is  possible  that  it  can  also  be  applied  to  larger  urban  areas  with 
greater  trip  end  opportunities  and  higher  travel  times.   Therefore, 
additional  verification  and  testing  should  be  conducted  on  a  few 
larger  urban  areas.   This  would  allow  for  testing  of  the  linearity  of 
origin  zone  amerage  travel  time  with  the  opportunity  origin  zone 
average  travel  time  at  higher  time  values. 

In  addition,  because  of  the  fact  that  one  lives  in  a  dynamic 
society  where  the  social  and  economic  environment  can  change  over  time, 
a  sampling  technique  should  be  developed  so  as  to  provide  a  method  for 
updating  or  validating  the  regression  equations  relating  origin  zone 
average  travel  time  with  the  opportunity  travel  time.   This  sampling 
technique  could  probably  be  performed  in  conjunction  with  the  limited 
household  trip  generation  survey  used  for  validating  or  updating  cross 
classification  trip  rates.   Recent  studies  have  shown  that  a  small 
sample  origin-destination  survey  can  probably  be  used  for  trip 
distribution  modelling  calibration  or  updating.   It  may  be  possible 
to  use  this  same  data  for  updating  the  relationship  of  origin  zone 
travel  time  and  opportunity  travel  time  which  is  used  in  the  proposed 
trip  distribution  procedure. 
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APPENDIX  A 
SPATIAL  LAND  USE  ACTIVITY  DISTRIBUTION  IN  LAFAYETTE, 
ANDERSON,  AND  MUNCIE,  INDIANA 

Spatial  land  use  activity  distribution  is  probably  the  greatest 
determinant  in  the  trip  making  decision  process.   In  the  gravity  model 
trip  distribution  step,  the  spatial  land  use  activity  distribution  is 
numerically  quantified  in  terms  of  the  network's  travel  time  matrix 
(skim  tree)  and  by  the  zonal  productions  and  attractions  by  trip 
purpose.   A  map  of  an  urban  area  along  with  a  general  description  can 
provide  additional  information  on  how  the  travel  pattern  may  be 
effected  by  the  land  use  distribution  and  various  topographical  features. 

Greater  Lafayette,  Indiana 

The  spatial  location  of  the  major  land  use  trip  generators  of  the 
Greater  Lafayette  area,  population  100,000,  in  combination  with  some 
significant  topographical  barriers  have  resulted  in  a  travel  pattern 
characterized  by  a  relatively  large  travel  time  variance  in  comparison 
to  other  small  Indiana  urban  areas.   A  map  showing  the  location  of  the 
major  land  use  trip  generators  and  the  topographical  barriers  of 
Greater  Lafayette  is  given  in  Figure  Al. 

Purdue  University,  located  on  the  west  side  of  the  Wabash  River, 
is  a  large  university  employing  approximately  9,000  persons  and  had  an 
enrollment  of  approximately  28,000  students  at  the  time  of  the  study. 
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Figure  Al .  Greater  Lafayette  Traasportation  and  Development  Study  Area 
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iL  is  the  largest  trip  generator  in  the  urban  area,  producing 
approximately  97,000  trips  per  day. 

The  city  of  West  Lafayette,  population  20,000,  also  located  west 
of  the  Wabash  River,  can  be  best  typified  as  a  bedroom  community  with 
no  major  industry  or  manufacturing  center.   The  central  business 
district  is  very  small  and  is  composed  of  small  service  oriented 
businesses.   A  small  shopping  area  is  located  at  tlie  northern  end  of 
the  city  and  another  along  the  Wabash  River.   The  city  of  West 
Lafayette,  in  itself,  is  not  a  major  trip  generator. 

The  city  of  Lafayette,  population  50,000,  located  on  the  east 
side  of  the  Wabash  River,  is  typical  of  many  older  cities,  consisting 
of  a  central  core  area  and  a  growing  ring  on  the  outskirts.   One  of 
Lafayette's  major  trip  generators  is  the  central  business  district 
(CBD) ,  located  near  the  river.   The  CBD  is  predominately  surrounded  by 
residential  areas.   The  residential  area,  in  turn,  is  bounded  by  a 
large  area  of  manufacturing,  wholesale,  retail,  and  commercial 
establishments  which  have  developed  along  the  U.S.    52  highway.   This 
area  is  still  growing  with  new  construction  and  scattered  subdivisions. 

Topographical  barriers  play  an  important  part  in  the  Greater 
Lafayette  area.   At  the  time  of  the  transportation  study,  only  four 
bridges  crossed  tlie  Wabash  River.   One  of  the  bridges,  due  to  old  age 
and  functional  obselescence ,  no  longer  exists.   The  importance  of  the 
river  as  a  topographical  barrier  is  best  exemplified  by  the  two  minute 
time  penalty  added  to  the  travel  times  on  each  bridge  to  correct  for 
overdistribution  in  the  gravity  model.   The  railroads,  also  shown  in 
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Figure  Al ,  are  also  topographical  barriers,  but  to  a  lesser  extent. 
No  time  penalties  have  been  assigned  to  links  crossing  the  tracks, 
however,  much  of  the  highway  network  layout  has  been  influenced  by  the 
location  of  the  railroad  tracks. 

This  scattering  of  the  three  major  trip  generators  and  the  topo- 
graphical barriers  therefore  play  an  important  role  in  the  determination 
of  Greater  Lafayette's  travel  pattern.   Unique  characteristics,  such  as 
these,  can  easily  complicate  the  transportation  modelling  process. 

Anderson,  Indiana 

A  map  of  the  Anderson,  Indiana  study  area.  Figure  A2 ,  reveals 
similarities  to  the  Greater  Lafayette  area.   Both  urban  areas  are 
divided  by  a  river  and  possess  a  university  or  college  within  their 
study  area.   Anderson's  study  area  population  of  approximately  90,000 
is  also  similar  to  Greater  Lafayette's  population.   At  this  point, 
however,  the  similarities  end.   As  was  found  in  the  calibration  of 
the  Anderson  gravitv  model,  travel  time  variance  is  approximately  one 
half  that  of  Lafayette,  Indiana. 

Anderson  College  is  much  smaller  when  compared  to  Purdue  University. 
Anderson  College  generates  approximately  8,000  trips  as  compared  to 
Purdue  University's  97,000  trips.   The  effect  on  travel  patterns  by 
Anderson  College  is  therefore  relatively  small. 

White  River  divides  the  Anderson  urban  area  into  a  northern  and 
southern  section.   As  compared  to  the  Greater  Lafayette  Area,  a  much 
greater  proportion  of  the  urbanized  area  is  located  on  one  side  of  the 
river,  therby  reducing  the  overall  effect  of  the  UTiite  River  as  a 
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topographical  barrier.   Secondly,  the  White  River  has  a  total  of  eight 
bridges  crossing  it  as  compared  to  Lafayette's  four  bridges,  thereby 
greatly  increasing  the  accessibility  to  and  from  the  northern  section. 

Another  factor  which  may  greatly  contribute  to  Anderson's  lower 
travel  time  variance  as  compared  to  Lafayette  is  the  fact  that  the  CBD 
and  two  large  manufacturing  centers  are  within  close  proximity  to  a 
centrally  located  high  density  residential  area.   This  is  in  contrast 
to  the  (Jreater  Lafayette  area  where  the  high  density  residential  areas 
and  major  trip  generators  are  spread  out  over  a  much  greater  area. 

Muncie,  Indiana 

Muncie,  Indiana,  as  shown  in  Figure  A3  with  a  study  area  popula- 
tion of  approximately  100,000,  is  also  divided  by  the  Wliite  River  and 
has  Ball  State  University  located  within  the  urbanized  area.   The 
Muncie,  Indiana  study  area  therefore  possesses  many  characteristics  of 
both  the  Greater  Lafayette  area  and  Anderson,  Indiana.   Muncie's 
travel  time  variance,  as  determined  from  gravity  model  calibration,  is 
approximately  equal  to  Anderson's  and  one  half  of  Lafayette's  travel 
time  variance. 

Ball  State  University  is  located  north  of  White  River.   Ball 
State  generates  approximately  34,000  trips,  excluding  intra-university 
trips,  and  is  therefore  a  major  trip  generator.   Another  indication  of 
the  importance  of  Ball  State  University  as  a  trip  generator  was  the 
necessity  to  use  K  factors  for  external-internal  trip  purpose  gravity 
model  in  the  original  transportation  study. 


126 


EDUCATIONAL,  BALL  STATE  UNIVERSITY 
CBD 

COMMERCIAL 

INDUSTRIAL,  TRANSPORTATION 


Figure  A3.   Delaware-Muncie  Metropolitan  Transportation  and  Urban 
Development  Study  Area 
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The  White  River  divides  Muncie  into  a  northern  and  southern 
section,  with  an  approximately  equal  proportion  of  land  use  activity 
on  each  side.   As  was  reported  in  Muncie's  transportation  study 
reports,  the  river  itself  is  not  a  major  topographical  barrier. 
Fourteen  bridges  cross  the  river  with  nine  of  these  being  in  the 
urbanized  portion  of  the  study  area.   Travel  time  penalties  were 
added  to  only  two  of  the  bridges.   This  situation  differs  greatly  from 
Lafayette's  and  may  therefore  help  to  explain  Muncie's  low  travel  time 
variance.   As  previously  noted,  Anderson's  bridge  situation  is  very 
similar  to  Muncie's  and  may  help  to  explain  the  similar  travel  time 
variances. 

A  review  of  Muncie's  land  use  map  shows  that  a  large  high  density 
residential  area  is  located  centrally  to  the  Muncie  CBD  and  within  a 
short  distance  to  Ball  State  University.   The  relatively  close 
proximity  of  these  three  major  trip  generators  result  in  trip  lengths 
of  a  relatively  uniform  length  and  therefore  probably  contribute 
greatly  to  Muncie's  low  travel  time  variance. 

Zone  and  Network  Comparison 
The  level  of  zone  and  network  detail  with  respect  to  population 
size  can  also  influence  model  results  and  the  statistical  measures 
used  to  test  model  accuracy.   The  three  urban  areas  of  Lafayette, 
Anderson,  and  Muncie  each  have  a  population  of  approximately  the  same 
size,  however  the  level  of  zone  and  network  detail  differ  significantly. 
Summarized  in  Table  Al  are  the  number  of  zones  and  one  way  links  for 
each  of  the  three  urban  areas.   An  analysis  of  the  effect  of  the  level 
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of  zone  and  network  detail  on  trip  distribution  and  traffic  assign- 
ment was  conducted  as  the  third  part  of  the  overall  study  on  synthetic 
travel  demand  modelling  (Hanscom,  E.  and  Sinha,  K.  C,  "The  Effects  of 
Simplifying  Traffic-Zone  and  Street-Network  Systems  on  the  Accuracy 
of  Traffic  Assignments  in  Small  Urban  Areas  in  Indiana",  Joint 
Highway  Research  Project  Report  No.  FHWA/lN/jHRP-79/15 ,  Purdue 
University,  1979. 


Table  Al.   Level  of  Zone  and  Network  Detail  Comparison  for  Lafayette, 
Anderson,  and  Muncle 


Internal 

External 

Total 

One-Way 

Urban  Area 

Zones 

Zones 

Zones 

Links 

Lafayette 

111 

17 

128 

1448 

Anderson 

169 

A3 

213 

2321 

Muncie 

271 

45 

216 

2667 
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APPENDIX  B 
CALIBRATION  OF  A  GRAVITY  MODEL  UTILIZING  REGRESSION  A1>IALYS1S 
AND  THE  ORIGIN-DESTINATION  SURVEY  TRAVEL  TIME  FREQUENCY 
DISTRIBUTION  INFORMATION 

Calibration  of  a  gravity  model  is  based  upon  the  repeated  adjust- 
ments of  a  set  of  friction  factors.   Adjustments  are  continued  until 
the  resulting  gravity  model  travel  time  frequency  distribution  satis- 
factorily approximate  to  the  information  obtained  from  the  home 
intervievv'  survey.   The  adjustment  technique  used  may  be  represented 
mathematically  as  follows. 

0D% 

where 

F  ,  =  friction  factor  value  at  time  t  for  iteration  k. 
tk 

0D%  =  percentage  of  all  trips  of  travel  time  t  as  obtained  from 

the  origin-destination  survey. 
GM%  =  percentage  of  all  trips  of  travel  time  t  resulting  from 

iteration  (k-1)  of  the  gravity  model, 
t    =  time,  usually  a  time  increment  of  one  minute. 
The  values  of  F  ,  are  then  plotted  against  time  t  and  either  a  smooth 
curve  is  drawn  to  represent  the  friction  factor  curve  or  regression 
analysis  is  used  to  fit  a  function  to  the  data  points.   This  appendix 
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discusses  the  use  of  regression  analysis  as  applied  throughout  this 
study. 

Calibration  of  the  gravity  models  for  Lafayette,  Indiana  and 
Anderson,  Indiana  were  perfoinned  using  a  negative  exponential  function 
to  represent  the  friction  factor  curve,  as  shown  in  Equation  B2. 

-3  t 
F  =  C  e  (B2) 

where 

F    =  friction  factor  at  time  t 

C    =  constant. 

3    =  calibration  parameter, 

t    =  time  Increment. 
The  transformed  linear  regression  equation  takes  the  form: 

)ln(F)  =  v.n(C)  -  Pt  (B3) 

where 

Hn  =   natural  log  function. 

F,C,3,t    =  as  previously  defined. 

Parameter  3  is  the  only  coefficient  of  importance  when  used 
In  the  gravity  model.   The  constant  "C"  changes  the  absolute  value 
only;  not  the  relative  value.   In  regression  analysis,  ?  would  be 
computed  using  the  following  equation: 

P  =ML^-I^  (B4) 

Yt^   -    {T.ty/rx 

where 

n   =  the  number  of  data  points. 

F,t  =  as  defined  previously. 
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Employing  Equation  B4  the  resulting  average  travel  times,  for 
both  the  Lafayette  and  Anderson  study  areas,  converged  to  a  value 
greater  than  that  derived  from  the  origin-destination  survey. 
Relatively  simple  statistical  reasoning  soon  remedied  this  problem. 
The  use  of  only  the  individual  data  points,  that  is,  one  F  value  for 
each  time  value,  infers  equal  occurrence  (importance)  to  each  data 
point.   In  simpler  terms,  it  implies  that  a  friction  factor  at  a  very 
low  time  value  or  at  a  very  high  time  value,  where  few  trips  occur, 
has  the  same  degree  of  importance  as  the  friction  factor  at  the 
average  travel  time  where  many  trips  occur.   Consequently,  the  data 
points  are  weighted  on  the  basis  of  their  frequency  of  occurrence. 
The  origin-destination  survey   travel  time  distribution  percentages 
are  thus  used  to  properly  weight  the  data  points.   The  appropriate 
equation  for  determining  the  value  of  the  Beta  parameter  can,  there- 
fore, be  written  as  follows: 

g  ^  Z(OD%  tF)  -  Z(OD%t)  ;:(OD%F)/100 
Z(OD%  t^)  -  (Z(OD%t))^/100 

where 

0D%  =  percentage  of  trips  that  occurred  with  a  trip  length  of 
time  t. 

Equation  B5  when  used  in  the  gravity  model  calibration  of 
Lafayette  and  Anderson,  Indiana,  produced  average  travel  times  satis- 
factorily matching  those  obtained  from  the  home  interview  surveys. 

The  FHWA  report  PLANPAC/BACKPAC  General  Information  makes  no 
mention  of  using  a  weighted  regression  analysis.   Recently,  however. 


132 


a  computer  program  called  AGM  (A  Gravity  Model)  has  been  added  to  the 
Urban  Transportation  Planning  System  (UTPS)  program  package  in  which 
a  type  of  weighted  regression  analysis  is  used  in  the  calibration  of 
gravity  models. 
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APPENDIX  C 
COMPUTER  PROGRAM  JMGM 

Computer  program  JMGM  was  written  for  the  purpose  of  using  the 
origin  zone  specific,  synthetically  self-calibrating  gravity  model  with 
the  LTMODEL  program  of  the  UTPS  package,  and  therefore  it  can  be  used  by 
most  planning  agencies.   A  PLANPAC  version  of  JMGM  was  also  coded.   Two 
other  computer  programs,  using  the  same  program  logic  of  JMGM  were 
also  written  for  usage  on  any  computer  system  with  a  FORTRAN  compiler. 
One  version  makes  extensive  use  of  core  storage  and  the  other  uses 
disk  storage  much  like  the  UMODEL  program.   Both  can  easily  be  modified 
for  data  input  and  output  conformity  and  therefore  can  be  integrated 
into  other  transportation  modelling  packages. 

Program  logic  flow  consists  of  five  major  steps  in  the  calibra- 
tion procedure: 

Step  1.   Distribution  of  trips  with  the  B  parameter  of  the 

negative  exponential  friction  factor  equation  equal  to 
zero  (all  friction  factors  equal  to  one)  and  subsequent 
computation  of  the  resulting  origin  zone  average  travel 
times. 
Step  2.   Distribution  of  trips  with  the  3  parameter  equal  to  a 

preset  value  and  subsequent  computation  of  the  resulting 
origin  zone  average  travel  times. 
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Step  3.   Using  the  origin  zone  average  travel  times  generated 
by  Step  1  and  the  developed  regression  equations, 
predict  the  origin  zone  average  travel  times  of  the 
calibrated  gravity  model. 
Step  4.   Using  the  data  points  from  Steps  1  through  3,  compute 

the  origin  zone  B  parameter  using  linear  interpolation. 
(Assuming  a  linear  relationship  exists  between  origin 
zone  average  travel  times  as  6  changes.) 
Step  5.   Distribute  trips  using  the  origin  zone  specific  B 

values  from  Step  4  and  compute  final  statistics. 
Figure  CI  shows  a  generalized  flow  chart  for  program  JMGM.   This 
flow  chart  along  with  the  program  listing  given  in  Figure  C2  with  its 
numerous  COMMENT  statements  provide  a  thorough  documentation  of  the 
program  logic.   Additional  information  can  be  found  in  Chapter  IV. 
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9  -  SIM  COmFUTED  AfJ^  =   SU'UIZ)  T(I,J) 

10  -  ZONAL  AUE.  at  BETA=0.0 

11  -  ZONAL  UhR.  at  EETA=0.0 

12  -  ZOiiAL  Ai.'E.  AT  BETA=PEETA(  TP :  OR  FINAL  UALUE 

13  -  ZONAL  UAR.  AT  BETA^PBETA^P )  OR  FINAL  'JALUE 

14  -  USED  O'iLY  AS  WORKING  STORAGE,   TEMPORARY  FRI 

INITIALIZE  ZONAL  MATRIX 

IF((ITNO.NE.n.OR.(IZ.NE.l))RETURN 
DC  3C  I=5rl3 
DO  30  J=1.0MAX 
Z(I,J)=0.0 


THIS  IS  THE  CODE  FOR  ENTRY  POINT  MDD13C 


REAL«4  PEETAC3).FCTR(39J 
PBETA(1)=0.15 
PEETA(2)=0,30 
PEETAi-3)=Q.30 


0* 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
t 
c 


THIS  IS  THE  CODE  FOR  ENTPY  POINT  MOB 1 3D 

ORIGiri  ZONE  SPECIFIC  S'vNTHETICALLY  SELF-CALIBRATING 

INTERNAL-INTERNAL  TRIPS  0NLY4. 


USER  PARAMETER  INPUT  'JIA 


=ARMS 


UPARMS(l)  =  1 
2 
3 

Ui=ARMS(2)  =  1 

2 


HBU 
HBO 
NHB 

COMBINED  REGRESSION  EQUATIONS 
LAFAYETTE  REGRESSION  EQUATIONS 


137=         JMGM 

1 

BY  MKE  FARRISJMGM 

2 

JMGM 

3 

JMGM 

4 

JMGM 

5 

JI'iGM 

B 

JMGM 

7 

JMGM 

8 

JMGM 

3 

JMGM 

10 

JMGM 

11 

JMGM 

12 

JMGM 

13 

JMGi-1 

14 

JMGM 

15 

JMGM 

IG 

JMGM 

17 

JMGM 

18 

JMGM 

13 

JMGM 

20 

JMGM 

21 

JMGM 

22 

CITON  FACTOR  JMGM 

23 

JMGM 

24 

J^GM 

25 

JMGM 

2G 

JMGM 

27 

JMGM 

23 

JMGM 

29 

JMGM 

30 

JMGM 

31 

»»»»»♦»»*»»*»»  j;-i;3f-| 

32 

»c»*->'*»*»»»«-»«jr'GM 

33 

JMGM 

34 

JhJM 

35 

JMGM 

3S 

JMGM 

37 

JMGM 

33 

JMGM 

33 

Jt'CM 

40 

JMGM 

41 

JMGM 

42 

JMGM 

43 

»»»»*»•«■»»»*-•*  JtlQH 

44 

♦«*♦-?-*»•»♦  K-»»»«j;-',Gf^ 

45 

JMGM 

46 

JMGM 

47 

JMCM 

43 

GRAUITY  MODEL  J^'GM 

49 

oMGM 

50 

JMGM 

51 

JMGM 

52 

JMGM 

53 

JMGM 

54 

JMGM 

55 

JMGM 

5S 

JMGM 

57 

JMGM 

JMGM 

59 

JMGM 

GO 

Figure  C2.   Program  Code,  JMGM,  for  Usage  in  UMODEL 
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C 
C 

c 
c 


c 
c 


101 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

r 


21 


100 
50 


S3 


3  ANDEPSOH  REGRESSION  EG'JATIQ'iS 

4  MUNCIE  REGRESSIOM  EQUflTIDNS 

r^QTE:   OMftX=DnAX=  NUMBER  OF  INTERNAL  ZONES  ocsunED. 
REf=lL»4  AL'E.  UfiR,  ATT,  TUAR.  SUM.  TRIP.  RAvE.  SX.  SXX,  TOT 
INITIALIZE  OR  CONTINUE 

IF((ITNO.EQ.l).AND.(IZ.Ea.l))GO  TO  101 

GO  TC  100 

COMPUTE  TOTAL  TRIPS 

TOT=:0.0 

DO  20  J=1.0MAX 

r0T=T0T+Z(2.J) 

INITIALIZE 

SX=0.0 

SXX=0.0 

ATT=0.0 

TUAR=0.0 

CIT£R=0 

AITEP=2 

DO  21  J=1.0MAX 

Z(8> J)=Zt3. J) 

Z(7.J)=T0T 

ITERATION  FLAGS  FOR  EECISION  TREE. 

CITEP=CALIBRATICN  FLAG 

0  -  EETA=0.0 

1  -  BETA=PBETA(TP) 

2  -  EETA=CrtLIERATED  UALUE  FROM  REGRESSION  EQUATIONS 

AND  INTERPOLATION. 

AITEP=AT-RACTICN  FLAG 

0  -  EEGirf  ATTRACTION  ITERATION 

1  -  CONTINUE  ATTRACTION  ITERATION 

2  -  ATTRACTION  ITERATION  COMPLETE.  I.E.  CLOSURE 

IFfIZ.EQ.i;uiRITE(S,50)ITNO.CITER.AITER 

FORNAT(10X.?HITNO   =    . 15, lOX, BriCITEP   =    , 15, lOX, 8HAITER   =    ,15) 

IFCtCITER.EQ.D.AND.i AITER.EQ.O).AND.(IZ.EQ.1))GO   TO   200 

IF((LlTEP.Ea.2),A>D.(AITER.En.0).AND.(IZ.Ea.l) iGO    TO   300 

IF(HlTER.EQ.0)GO   TO   400 

IFCAITER.EQ.IJGO  TO  500 

GO  TO  GOO 


C 

C. 

C. 

C 

C 

C 


CITER=1,AITER=0,IZ=1    SET  BETA(I)*S  AND  FRICTION  FACTORS  TO 

200  TP=UPARMSi'l) 
DO  201  J=1,2S 

201  FCTR(J)=EaP(-P3ETA(TP)»J) 
GO  TO  SS 


JMGM 
J'^GM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
Jf^CM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 

j.-;gm 

JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 
jMGM 
JMGM 

jr-.GM 

JMGM 
JMGM 
JMGM 
JMGM 

JMGM 

^r.GM 

JMGM 
GIUENJMGM 
JMGM 
JMGM 
JMGM 
JMGM 
JMGM 


61 

62 

63 

64 

65 

66 

67 

68 

63 

70 

71 

72 

73 

74 

75 

76 

77 

73 

79 

80 

81 

82 

S3 

£4 

£5 

26 

87 

83 

ES 

SO 

91 

32 

S3 

S4 

55 

9S 

57 

S3 

93 

10  0 

101 

102 

103 

104 

105 

106 

107 

103 

109 

110 

HI 

112 

113 

114 

115 

116 

117 

lis 

113 
120 


Figure  C2.  Continued 
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C 

c. 

c. 

c 

r 
C 


CITES=2.fiITER=0, IZ=1 


FIND  BETACI.'  S  FROM  REGRESSIGN  EOUA. 


C 
C 
C 
C 
C 

c 
c 
c 
c 
c 

c 

r 

c 

c 


c 

c 


c 
c 

G 
C 

c 

c 


L 

c 
c 


300  RTP=U?ARnS ( 1 ; + ( UPRRMS ( 2 ) - 1 ) *3 
ATT=0.0 
TUPf?=O.C 
SX=0.0 
5Xo=0.0 

DO  IS  j=.'.,onft.x 

IF(ZC2.  J^.EQ.O.OGQ  TO  13 

COMPUTE  nAUE.   AUERftGE  FROM  REGRESS  I  ON  EQ'JATIDM 

AUE   =  ZONAL  AiJERAGE  TRAUEL  Tir-^E. 

AUEO  =  ZO;'IAL  h'.'ERAGE  TRAUEL  TI'-IE  AT  BETA  =  0.0 

UAR   =  ZOHAL  T;?AUEL  TIME  UARIh;-ICE  AT  CALlERATEIi  BETA  UALUE. 

UARO  =  ZONAL  TRAUEL  "^IHE  UARIAHCE  AT  BETA  =  0.0 

TERM  =  AUERAGE  ORIGIN  ZONE  TERMINAL  TIME. 


GO  TG  (1.2,  3,  4, 5.  B,  7. 8.  9, 10.  11.  12),  RT? 


COMBINED  REGRES3IGN  EGLJfiTIONS-LAFAVETTE.  ANDE=30N.  AND  NLNCIE  IND. 

HBU   R»*2  =  0.SE5 

RAUE=-0.4327-^1.2rSl»Z(4,  J)H-0.3921'*(Z(10.  J)-Z(4,  J)) 

GO  TG  19 

HBO   R»*2  =  0.877 

RAUE=-0.7433-^1.3543*Z(4.J)+0.7557*(Z(10.Ji-Z(4.J)) 
GO  TO  13 

NHB   R**2  =0.735 

RAUE=0.5245+1.341E»Z(4.J)+0.S073*(Zri0.J)-ZC4,J)) 
GO  TO  13 


LAFAYETTE,  INDIANA. 

HBUl   R**2  =  .893    AREAUIDE  UfiR=22.99;   UAP0=24.99 
AUE=10.60:   AUE0=12.32;   VAR=12.17:   UARC=14.00;   TERM=1.08 
RAUE=-l,0331  +  l.C3343Z(4,J)-r0.3328»(Z(10.J)-ZC4,  J)) 
GO  TO  19 


KEO   R»*-2  =  .853    AREAUIDE  UAR=12.G3:   UAR0=25.5-i 
AUE^8.55;   AUE0  =  11.78:   UriR=10.78;   UAR0=15.B3;   TERr 
RAUE=-l.b7Sl  +  1.145E*Z(4,  J)+0.S39r»tZ(10.  JJ-Z';4.  jn 
GO  TO  19 


=  1.11 


NHB   R»*2  =  .759    AREAUIDE  UAR=14.03;   UAR0=17.75 
AUE=9.12:   AUE0  =  11.42:   l'AR=3,34;   UAR0=12.91:   TER''=1.23 
RAUE=0.2131+l.c323*Z(4, J>+0.S744>(ZC10. J)-Z(4,J)) 
GO  TO  19 


ANDERSON,  INDIANA 


JMGM 

121 

.  J.rGM 

122 

,JMGM 

123 

J:-'GM 

124 

J;1GM 

125 

JMGM 

12S 

JMGM 

127 

JMGM 

123 

JMGM 

129 

JMGM 

130 

JNGM 

131 

JMGM 

132 

JMGM 

133 

JMGM 

154 

JMGM 

135 

JMGM 

13S 

JMGM 

137 

JMGM 

133 

JMGM 

133 

JMGM 

140 

JMGM 

141 

^."IGM 

l'i2 

JMGM 

1-3 

J-.GM 

144 

JMGM 

I'^S 

JMGM 

14^ 

,  JM^M 

147 

'jMGM 

143 

jMGM 

143 

uMGM 

15C 

JMGM 

151 

jMGM 

152 

JMGM 

153 

JMGM 

154 

JMGM 

155 

JMGM 

155 

Jf^GM 

157 

JMGM 

158 

JMGM 

159 

JMGM 

IGO 

JMGM 

Icl 

JMGM 

1P2 

JMGM 

153 

JMGM 

1E4 

JMGM 

165 

JMGM 

lEG 

JMGM 

1E7 

J.'iGM 

1G8 

J^'GM 

1E3 

JMGM 

170 

jMGM 

171 

JMGM 

172 

JMGM 

173 

JMGM 

174 

JMGM 

175 

JMGM 

172 

JMGM 

177 

JMGM 

178 

JMGM 

179 

JMGM 

130 
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C 
C 

c 


c 
c 


c 
c 
c 


c 
c 
c 


HEU        R»-»2   =    -873  AREPUIDE  UAR=9.B1;      UAR0=9.51 

fiUE=11.43;      ftiJE0=13.10;      UnR=S.53;      U0=7.12:       IERri=2.38 
R^UE=-0."650  +  1.130£^2(4. J)  +  0.94G?'*(Z(10.  J'-Z(4.  J)) 
GO   TO    13 

HBO        R»*2   =    .RS4  fiREfi^IEE   U^R=S.5l:       UPiRO^ll.lO 

aUE=10.C4:      fHi;EC  =  lE.30:      'JPlR=7.a4:      UAR0=S.20:      TER"=2.37 
RAUE=1.0Ell+i.0G40»Z(4, J)+O.S504»fZaO. J)-ZC4, J)) 
GO   TO    IS 

NHB        R»-^2=    .210  APEAUIEE   UAR=7.S2:       UAR0=9.20 

fiUE-=9.50:       P,IE0=11.7P;       UPR=R.35;       iJPR0=7.i2;       TERr!=2.B4 
RAl'E=27455B+1.0SG3»Z(4.J)--0.4£80*CZ(10.  J)-Z(4.  J)  ) 
GO   ^0    19 


ML^'CIE,    INDIAf-'A 

riEl-J        R»*2   =    .772  AREAl-OIEE   UAR=i:.lS:       UAR0=11.51 

Pl'E=9.07;       AUE0=9.S2;       Uft"=7. 03;       UP!R0=7.1i;       TERri=i.33 

10  RAUE=0.8E1S+0.7^43»Z':4.  J")+0.8491»(Z(10.J)-Zi'4.  J)) 
GO   TO    19 

HBO   R«*2  =  .743    AREfiUIEE  UCiR=9.E7;   ';fiR0=11.83 
ril.'E=7,48;   AU£0=S.53:   UAR=S.40;   iJfiR0=7.79;   TERM=i.34 

11  PAUE=0.1£;03+1.0595--2(4,  Ji+0.71G2-*(Z'10.  J)-Z':4.  J)  > 
GO  TO  19 

riHB   R*»2  =  .505    AREAUIDE  UAR=7.S9;   UAR0=9.23 
mUE=G.93;   AfE0=8.25;   UftR=5.G5;   Uf!R0=6,B5;   TERM=1.49 

12  RAUZ=0.7980+1.4132*Z(4. J)+0.53d4»(Z'10,J)-2(4.J)) 

19  CONTINUE 

CCMPUTE  BETA.  LIhEPR  IMTE-^OLATIDN. 

2(G.J)  =  (PEETfi(TP)/'(Z(:i2,  Ji-ZriO.  J)))*(Rfi'OE-Z(10,J)) 
IF(Z(12.J).GT,Z(10. J))Z(S, J1=0.0 
IFCZCG. J).GT.0.5;Z(G.J)=0.5 

IS  COMTIMUE 
GO  TO  33 


P.ITER=0.IZ=1 

400  DO  401  J=1.DMAX 
Z(8. J)=Z(3.Ji 

401  Z(9-J)=U.O 
AITEP=1 


INITIALIZE  BEFORE  ATTRACTION  ITERATION  BEGINS. 


C 
C 
C 

c 
c 


BEGIN  ATTRACTION  ITERATION. 
SUMJ  R'J.1*F(I,J) 
500  Z(7.-IZj=C,0 

DO  501  j=i,ona;< 

IF(IZ.E0.jnNT(J,n=Z(5.J) 


JMGM 

1?1 

jncr 

132 

jnor 

153 

j'iGn 

1S4 

jriGM 

lo5 

jr.cr^, 

18G 

JMGM 

187 

JMGM 

188 

JMGM 

139 

j;iGM 

190 

jrcM 

131 

JMGM 

192 

Ji^GM 

1S3 

JfrGM 

194 

JMGn 

195 

jr.GM 

ISG 

^■;iGi1 

137 

jriGN 

133 

J^iGM 

153 

jnon 

200 

jriGM 

201 

J'-iGM 

202 

J1GM 

203 

JMGM 

204 

JMGM 

205 

jroM 

20G 

jrGM 

207 

JiTGN 

203 

JMGM 

203 

JMG1 

210 

jt^Gn 

211 

jf-iGr 

212 

JMGM 

213 

JMGM 

214 

jr-'GM 

215 

JNGN 

21£ 

JMGll 

217 

JMGM 

218 

J.IGM 

219 

J-^Gfl 

220 

j^.zr, 

221 

JMGM 

222 

J^GM 

223 

JMGM 

224 

JMGM 

225 

JMGM 

225 

Jf^GM 

227 

223 

JMGM 

223 

JMGM 

230 

JMGM 

231 

JMGM 

232 

JMGM 

233 

JMGM 

234 

J'-GM 

235 

JMGM 

23G 

JMGM 

237 
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233 
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240 
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IFdNTCJ.  D.EO.O.OJGO  TO  501 

Za4,  J)  =  1.0 

ADD  TERriliNAL  TIMES  TO  SKIM  TREE  TRWEL  TIME. 

TiME  =  INT':j.l)-Z(4.  IZ)+Z(4,J) 

IF(T:,1E.GT.39)TIME=SS 

IF(CITER.£a,ljZ(14. Ji=FCTP(TIME) 

IF(CITER.E0.2)Z(14,Jj=EXF(-2(S,IZ)*TIt1E) 

SAF  =  SAF  +  A  •  F 

ZC7.  iz)=zcr,iz)+z(8.j)»zn4.j) 

501  ccriTir<uE 

SL'M   COMPUTED   A-:-TRACTIONS 

cc  502  j=i.oma;< 

IF(IMT(J,1).EQ.O.O)GO   TO   502 

51^1   =   5P)   +   p   »   rt   »   F   /   SAP 

ZO.  J)=Z(3.  J)+2(2,IZ)*2(8,J)»Z(14,J)/Z(7,  IZ) 

502  COMTI.-iUE 


C 

C 
C 

c 


IF  IZ.EC.OMAX. 
IF  IZ.LT.GMAX. 


TEST  FOR  CLCSL'^E 
CCMTIN'JE  ITERft~IOM 


C 

c 
c 

c 
c 
c 


IF(IZ.EQ.OMAX)GO  TO  ?00 
RETURh  1 

COMPUTE  ERROR  FOR  CLOSURE  "EST 

?00  SUM=0.0 

CO  ?'01  J=1,CMAX 
701  3UM=SUM+ABS(Z(3,J)-Z(3,J)) 

SUM=SUM/TCT 

OME  PERCEMT  CLOSURE  TEST 

IFCSUM.LT.O.ODGO  TO  702 

CLOSURE  TEST  FAILED,   COMPUTE  MEU  ATTRfiCTIOISS. 

DC  703  J=1,0MAX 
IF(Z(8. J).EQ.0.0)GO  TO  703 
Z(8. J)=Z(3, J)»Z(8. J)/ZCS- J) 
703  Z':9.J)=0.0 
RETURN  1 


CLOSURE  TEST  PASSED.   SET  AITER=2  FOR  AUE.  AMD  UAR.  ITERATION. 

702  AITER=2 

IF  C:TER=2,  SET  ITr-a=ITER-l  FOR  LAST  AUE.  UAR.  ITERATION. 

IF(CITER.EQ.2)ITr'0=ITER-l 
RETUF1  1 


AUERmGE   AMD   UAPIAI^CE   COMPUTATION 
INITIALIZE   ZONAL   AUE.    AND   UAR. 
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GOO 

ftUE=0.0 

JMGM  301 

UAR=0.0 

JMGN  302 

DC  533  J=1.Z0ME3 

JMGM  203 

553 

TROUT(J,1)=0 

JMGM  304 

IF(((CITER.EQ.0).OR.(CITER.EQ.2)).flND.(lZ.EQ.l))Ga  TO  538 

JMGM  305 

GO  TO  537 

JMGM  30e 

533 

fiTT=C.G 

JMGM  307 

TUP3=0.0 

JMGM  303 

SX=0.0 

JMGM  303 

SXX^O.O 

JMGM  310 

c 

JMGM  311 

c 

BEGIN  ^UE.  AND  UAR.  COrPUTPlTIO'H 

JMGM  212 

c 

JMGM  313 

c 

IF  GPIGIM  ZDME  PRODUCTION=0.  SKIP  fi'JE.  AND  UAR. 

JMGM  314 

c 

JMGM  315 

537 

IF(~:2.I2).LE.1.0^GO  TO  G05 

JMGM  31S 

DO  SCI  J=1,GMAX 

JMGM  317 

c 

JMGM  313 

c 

COM'^UTE  CELL  IHTEPCHAMGE,  I.E.  TRIPS. 

JMGM  313 

c 

JMGM  320 

IF(IZ.Ea.J)IMT(J.l)=ZC5,J) 

JMGM  321 

IFCIMTCJ. l).Ea.0.0)GO  TO  GOl 

JMGM  322 

Z(14. J)=l 

JMGM  323 

TIME=IMTrj, 1)+ZC4.IZ)+Z(4.  J) 

JMGM  324 

IF(TIt1E.GT.39jTIME=93 

JMGM  325 

IF(C!TER.E0.1)Z(14, J)^FCTR(TIPE) 

JMGM  32S 

IF ( CI TER.  EQ.  2 ) Z(  14.  j)=EXr-( -Z(G,IZ)»TIME) 

JMGM  327 

c 

TRIP  =  P  »  A  •  F  /  SAF 

JMGM  323 

TRIP=Z(2.  IZ)*Z':S.  J)*Z'14.  J)/Z(7,  IZ) 

JMGM  329 

c 

JMGM  330 

c 

SAUE  TRIP  IN  TROL'T  IF  LAST  ITERATICM.  I.E.  CITER=2 

JMGM  331 

c 

JMGM  332 

IF( (CITEP.E0.2) .AND. (AITER.EQ.2) ) TROUT (J,  1)=TRIP 

JMGM  333 

c 

JMGM  334 

c 

SUN  ZONAL  T,?AUEL  TIMES. 

JMGM  335 

c 

JMGM  35G 

AUE=fiUE+TRIF»TINE 

JMGM  337 

SOI 

CONTINUE 

JMGM  333 

c 

JMGM  333 

c 

ADD  ZONAL  CQNTRIEUTIOli  TO  AREAl^IDE  AUERAGE  TRAUEL  TlflE 

JMGM  340 

c 

IF((CITER.EO.0).OR.CCITER.EQ.2)) 

JMGM  341 

r 

JMGM  S-iE 

IF((CITEP.EQ.O).OP.CCITER.EQ.2  5)ATT=ATT+AUE/TOT 

JMGM  343 

c 

JMGM  3^4 

c 

COMPUTE  ZONAL  AUERAGE  TRAUEL  TINE 

JMGM  345 

c 

JMGM  34S 

ftUE=AUE/Z(2.IZ) 
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